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ABSTRACT 

Galaxies are not limited to simple spheroid or bulge -i- disc morphologies. We explore the 
diversity of internal galaxy structures in the Coma cluster across a wide range of luminosi¬ 
ties (—17 > Mg > —22) and cluster-centric radii (0 < rduster <1-3 r 20 o) through analysis of 
deep Canada-France-Hawaii Telescope i band imaging. We present 2D multi-component de¬ 
composition via GALFiT, encompassing a wide range of candidate model morphologies with 
up to three photometric components. Particular focus is placed on early-type galaxies with 
outer discs (i.e. SOs), and deviations from simple (‘unbroken’) exponential discs. Rigorous 
filtering ensures that each model component provides a statistically significant improvement 
to the goodness-of-fit. The majority of Coma cluster members in our sample (478 of 631) 
are reliably fitted by symmetric structural models. Of these, 134 (28%) are single Sersic ob¬ 
jects, 143 (30%) are well-described by 2 component structures, while 201 (42%) require more 
complex models. Multi-component Sersic galaxies resemble compact psuedobulges (n ~ 2, 
Re ^ 4 kpc) surrounded by extended Gaussian-like outer structures (Re > 10 kpc). 11% 
of galaxies (N = 52) feature a break in their outer profiles, indicating ‘truncated’ or ‘anti¬ 
truncated’ discs. Beyond the break radius, truncated galaxies are structurally consistent with 
exponential discs, disfavouring physical truncation as their formation mechanism. Bulge lu¬ 
minosity in anti-truncated galaxies correlates strongly with galaxy luminosity, indicating a 
bulge-enhancing origin for these systems. Both types of broken disc are found overwhelm¬ 
ingly (> 70%) in ‘barred’ galaxies, despite a low measured bar fraction for Coma (20 ± 2%). 
Thus, galaxy bars play an important role in formation of broken disc structures. No strong 
variation in galaxy structure is detected with projected cluster-centric radius. 

Key words: galaxies: clusters: Abell 1656, galaxies: elliptical and lenticular, cD; galaxies: 
evolution; galaxies: formation; galaxies: structure 


1 INTRODUCTION 

Lenticular (SO) galaxies occupy the crux of the Hubble se¬ 
quence, representing the morphological intermediate between disc- 
dominated spiral galaxies and spheroidal ellipticals. However, it 
remains unclear whether SOs are evolutionary intermediates be¬ 
tween star-forming late-type galaxies and mainly passive early-type 
galaxies (ETGs). This evolutionary link has been extensively in¬ 
vestigated with emphasis on the transformation of spirals into SOs 
via quenching their star formation (see |Barr et al.|2007| |Arag6n-| 
|Salamanca|2008[|Barway et al.|200^ . 

Classically, SOs comprise a spheroid-shaped bulge component 
and a smooth disc with little or no interstellar dust or star formation. 

JL 

These bulge and disc structures are well described by Sersic (r " ; 
|Sersic|196^ and exponential profiles respectively. Conversely, gi¬ 
ant elliptical galaxies are traditionally viewed as smooth, single 
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spheroid systems well-described by a de Vauccouleur’s profile 
(Sersic n = 4; |de Vaucouleurs|1948^ . The morphological distinc¬ 
tion between these two classes can be unreliable depending on disc 


Strength, galaxy inclination, or observation depth (|Kent 

1985||Rix| 

|& White] 199011 Jprgensen & Eranx|1994[|van den Bergh 

2009bV 


[van den Bergh| ( [7976l l introduced the idea that the SO morphol¬ 
ogy encompasses multiple distinct classes of galaxy (SOa-c; analo¬ 
gous to the spiral Sa-c types), differing in luminosity and evolution¬ 
ary pathway (see also |van den Bergh|1990||2Q09a^. This concept 
was supported by kinematic studies of SOs (e.g. Dressier & Gunn| 
|1983| l, which demonstrated equivalence of the rotational properties 
of disks in SOs and spiral galaxies. More recently, t his idea has been 
developed further by the ATLAS®° group (e.g. Emsellem et al. 


|201 l[|&ppellari et al.|20l7] see also |Kormendy & Bender|2012^ . 
In this paradigm, most ETGs form a continuous sequence of rotat¬ 
ing, quiescent galaxies with specific angular momentum increasing 
with Hubble T stage. Hence, with sufficient signal-to-noise (S/N), 


© 2015 RAS 





























2 J.T.C.G. Head et al. 


discs should be detectable in many galaxies classically typed as el¬ 
liptical. 

With increasing local environment density, the morphologi¬ 
cal fraction of galaxies becomes increasingly dominated by ETGs 
(particularly SOs). Conversely, spiral galaxies are rare in the dense 
cluster environment. This morphology-density relation ( |Dressler| 
|1980[ l implies that the cluster environment plays an important role 
in the evolution of SOs from spirals (or spiral-like progenitors). The 
mechanisms potentially responsible for this evolution (see review 
in |Boselli & Gavazzi||2006l ) can be broadly categorised as disc¬ 
fading (e.g. gas-stripping) or bulge-enhancing (e.g. tidal interac¬ 
tions/mergers). While the latter category is traditionally thought 
of as disc-disruptive, it has been demonstrated that SO morpholo¬ 
gies can survive merger-based quenching ^Eliche-Moral et al.|2013[ 
IQuerejeta et al.|20T4) . 

The well-studied Coma cluster (Abell 1656) possesses one of 
the richest ETG populations in the local universe. As such. Coma 


As noted in |Dullo & Graham| ( |2014^ , care must be taken to report 
“real” structural components, rather than overfitting galaxies with 
unnecessarily complex models. 

Here, we build upon these pioneering studies by characterising 
the multi-component internal structures of galaxies within a wide 
radial area (0 < Tciuster < 1.3r2oo) of the Coma cluster (and an 
absolute magnitude range—17 > Mg > —22) using deep Canada- 
France Hawaii Telescope (CFHT) i band imaging data. 

The decomposition analysis reported in Paper I is extended 
by using a wider suite of candidate models (including 2- and 3- 
component broken disc galaxies) in order to explore the diversity 
of galaxy shucture in the Coma cluster. Thus, we reinvestigate the 
structural morphologies of all Coma cluster galaxies investigated in 
Paper I, including the ~ | previously removed from analysis as not 
well-described by an archetypal bulge -l- disc model. While a pri¬ 
mary goal of this analysis is the investigation of Type I, II, and III 
discs galaxy structures, the extended range of (multi-component) 
models is necessary to avoid mis-classihcation of additional com- 


1902|[Shapley|1934|[Dressler|1980| van Dokkum et al.|2014| and ponent structures (e.g. bars or rings) as surface brightness profile 

characteristics of ETGs (e.g.lLucey et al.|1991| Bower et al.|1992 

breaks. Bayesian model selection and sample filtering are applied 

Jprgensen 

19991|Hudson et al.|20101|Lansbury et al.|2014 Weinzir 

] to avoid overfitting, and to ensure that best fit models are reliable 

et al.|2014 

J. In addition. Coma encompasses a wide range of envi 

representations of the underlying galaxy structures. 


ronment conditions (~ 100 x difference in galaxy density between 
the core and the virial radius), allowing in-depth investigation of 


radial trends of environment-mediated processes l Gavazzi| | _ 

[Guzman et al.|1992[[Carter et al.|2008[|Gavazzi et k.|2010||Smith| 

[et al.|2012|[Cappellari|20T3l[Rawle et al.|2013[ l. 


In 


Head et al. 


(|2014[ hereafter ‘Paper P), we presented bulge 


- disc decompositions of ~ 600 Coma cluster galaxies, demonstrat¬ 
ing that ~ I of Coma ETGs are well-described by an archetypal 
SO (central Sersic bulge -l- outer exponential disc) model morphol¬ 
ogy. Focusing exclusively on these archetypal galaxies, we found 
that bulges of SO galaxies resemble pseudobulges (n ~ 2, i7e ~ 1 
kpc), while their discs were measured to be intrinsically smaller, or 
brighter than equivalent structures in star-forming spirals. A bulge 

— disc colour separation of ~0.1 mag was measured in g — i 
(~0.2 mag in u — g), indicating either a ~2-3x age difference, 
or a ~ 2x metallicity difference between these components. Nev¬ 
ertheless, both components were found to contribute to the galaxy 
red sequence (colour-magnitude) trend. 

Evolutionary pathways will not necessarily preserve the 
archetypal SO morphology. Furthermore, the simple exponential 
model (Type I; [Freeman] [1970^ adopted in most decomposition 
studies does not fully represent the observed range of SO outer 
disc structures. ‘Broken’ disks have been observed for SO and spiral 
galaxies ( |Freeman|1970[|Erwin et al.|2008^ , wherein surface bright¬ 
ness profiles beyond a break radius deviate either downwards (i.e. 
fainter; ‘Type IF) and upwards (i.e. brighter; ‘Type IIP) relative to 
a simple exponential (‘Type P) profile. Such profiles result from the 
redistribution of stars due to evolutionary processes. For example, 
truncated discs may be formed when stars are physically removed 
from a galaxy’s outer regions (e.g. during tidal interaction), while 
anti-truncated discs may result from merger events ( [Youn g er et al.[ 
[2007[ [Borlaff et al.[[2014) . Thus, investigation of galaxies with a 
wider range of structural morphologies provides a more complete 
picture of the ETG formation mechanisms. 

Previous investigations of multi-component ETG structures 


(e.g. Michard|1985 Capaccioli et al. 

1991| Laurikainen et al. 

2005 

Janz et al.|2012||Huang et al.|2013 . 

anz et al.|2014|[Weinzir 

et al. 


2012 


We investigate four main questions regarding galaxy evolu¬ 
tion: Does the multi-component structure of giant ellipticals sug¬ 
gest the ‘puffing-up’ of a compact progenitor, or the accumulation 
of additional structures around a compact spheroid? Are broken 
discs structures (truncated or anti-truncated) correlated with the 
properties of the bulge/bar components? Do the structures of Free¬ 
man Type II galaxies indicate physical truncation of discs? Does 
such a truncation scenario explain the apparent size offset of SO 
discs relative to star-forming spirals reported in Paper I? 

The structure of this paper is as follows: first, in Section 2 we 
summarise the MegaCam imaging data and galaxy sample selection 
criteria used in this work. Secondly, Section 3 describes the multi- 
component decomposition methodology, highlighting differences 
from the bulge-disc decomposition pipeline previously-described 
in Paper I. Thirdly, in Section 4 we present the resulting galaxy 
morphology (model) fractions, including a census of disc types. 
Furthermore, we explore the properties of galaxies comprising mul¬ 
tiple distinct Sersic structures, and galaxies containing disk breaks. 
Finally, a discussion of possible formation pathways for broken 
discs is presented in Section 5. 

Throughout this paper, we make use of the following notation 
conventions: Fitted model structures (see Section]^ are indicated in 
italics (e.g. ‘S’ for pure Sersic) to distinguish them from morpho¬ 
logical classifications (e.g. ‘SO’). Disk break types (i.e. Freeman 
types; untruncated, truncated, anti-truncated) are denoted with Ro¬ 
man numerals (e.g. ‘Type IP), and galaxies containing such struc¬ 
tures are referred to as Type I, Type II, or Type III galaxies. Con¬ 
versely, galaxy types using Arabic numerals (e.g. ‘Type 2’) refer to 
[Allen et ‘ak] ( [2006} surface brightness profile types (see also Sec- 
tion[^. The Type 1 profile is a special case describing a central 
bulge and an outer (exponential or broken exponential) disk, and 
is referred to as an ‘archetypal SO’ profile (‘archetypal’, or ‘SO’ as 
shorthand). All other [Allen et al.[ j2006^ types are referred to as 
‘atypical SO’ profiles (or simply ‘atypical’). 

We use the WMAP7 cosmology: Hq = 70.4kms“^Mpc“^ 
(i.e. /i 7 o = 1.01), Dm = 0.272 and Da = 0.728 ( [Komatsu et al.[ 


2011 1 . Using 2 cMB(Coma) = 0.024, the luminosity distance for 


[Laine et al.|2014|l are typically limited by (relatively) small 


the Coma cluster is 104.1 Mpc, and the distance modulus, m — 
M = 35.09. At this distance, P corresponds to 28.9 kpc. Taking a 


galaxy samples from narrow fields of view or ID profile analyses. 


value for velocity dispersion of acoma = 1008 km s (jStruble & 
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Rood 


et al. 


1999' and virial mass, M 200 = 5.1 x Gavazzi 


20091, the virial radius, r 2 oo, for Coma is 2.2 Mpc (~ 75'). 


2 DATA AND INITIAL SAMPLE 

This study makes use of the data as previously described in Pa¬ 
per I. To recap: optical imaging covering a total of 9 deg^ of the 
Coma cluster in the i band was acquired using the MegaCam in¬ 
strument on the 3.6 m CFHT during March - June 2008 (Run ID 
2008AC24, PI: M. Hudson). Total (coadded) exposure times of 300 
s were obtained for each observed field, yielding ~ 12 x deeper 
imaging data (from D^fexp) compared to SDSS (2.5 m telescope, 
53 s exposures). The MegaCam frames were sky-subtracted during 
pre-processing using a 64 pixel mesh. A point spread function (psf) 
full-width half-maximum (fwhm) of between 0.65” and 0.84" was 
typical. The pixel scale was ~ 0.186 arcseconds/pixel. 

The initial sample for analysis was selected from SDSS (DR9) 
catalogue galaxies in the 3 deg x 3 deg (= 5.2 Mpc x5.2 Mpc, 2 
r 2 oo X 2 r 2 oo) area covered by the MegaCam observations. A limit 
of —17.1 > iVfJ^was applied to ensure sufficient signal-to-noise 
(S/N) for reliable measurement of galaxy bulge and disc structures. 
These targets were limited to the redshift range 0.015 < z < 0.032 
(heliocentric vcoma ± 2.5(Jid) to ensure that only cluster members 
were included. Unlike Paper I, no colour cut is made during sample 
selection. Thus, an additional ~60 blue galaxies (g — r ^ 0.5) 
are included in the present work, yielding an initial sample of 631 
Coma cluster members. 

To illustrate the initial image analysis undertaken, and the di¬ 
versity of photometric structures observed, we show in Figure[T]the 
major axis surface brightness profiles for six representative Coma 
galaxies. As well as displaying the profile derived from the stacked 
image data (black points), we show the results from the individ¬ 
ual MegaCam exposures (grey lines). The radial limit used in our 
profile analysis is shown (grey vertical bars)', this corresponds to 
1 — 3% of the sky brightness level (i band ~ 25 mag arcsec“^). 
Within the fitted area, there is very good agreement between the 
stack and individual image profiles. This demonstrates that the 
stacking process, including the choice of sky-grid mesh size, has 
little or no effect on the derived profiles and hence on the 2D sur¬ 
face fitting analysis reported below. 

Here we briefly note some of the key features apparent 
in the surface brightness profiles for these six galaxies. Galaxy 
1237667322723369088 has a distinct break in profile at ~ 8" 
with a second, downward-bending exponential-like outer shape. 
1237667444048593359 has a slight upward-bending outer struc¬ 
ture, although this may be affected by the nearby contaminating 
galaxy. 1237667323797504020 displays a relatively weak interme¬ 
diate exponential-like structure, and an outer downward-bending 
shape. 1237667444048527399 has a very distinctive upward- 
bending, exponential-like outer structure. 1237667444585595093 
has a weak, downward-bending exponential-like outer structure. 
1237667442974392369 has a weak, upward-bending outer struc¬ 
ture. 

In Section 4 the model surface brightness profiles for these 
six galaxies derived from the multi-component fits are shown over¬ 
plotted on the data points. Work in progress will provide a de¬ 
tailed comparison of the surface brightness profiles of Coma clus- 


^ Note that while the sample is defined based on g band photometry, we 
analyse the i band data in this paper. 


Component 

Parameter 

Description 

D 

rtii 

Total i band magnitude 


Rs 

Exponential scale length 


q 

Axis ratio {b/a) 


4> 

Position angle 

BorS 

m 

Total i band magnitude 


Re 

Effective half-light radius 


n 

Sersic index 


q 

Axis ratio (b/a) 


$ 

Position angle 

C 

m 

Total i band magnitude 


Re 

Effective half-light radius 


n 

Sersic index 


q 

Axis ratio (b/a) 


$ 

Position angle 


Co 

Boxiness/Diskiness parameter 

Dd 

Mbrk 

Surface brightness at r^rk 


-^s,in 

Inner exponential scale length 


Rs,out 

Outer exponential scale length 


rbrk 

Break radius 


q 

Axis ratio (b/a) 


$ 

Position angle 


Table 1. Table of the model components used during decomposition anal¬ 
ysis, including descriptions of their free parameters. Note that the ‘bulge’ 
label (B) is used to describe the central Sersic component. All components 
in a model share centroid position parameters (x, y). Cq is defined in|Peng| 

Fnripoio) . 


Model 

Label 

k 

^comp 

Sersic 

s 

1 

1 

Sersic + exponential 

BD 

11 

2 

Boxy Sersic + exponential 

CD 

12 

2 

Double Sersic 

BS 

12 

2 

Sersic + broken exponential 

BDd 

13 

2 

Boxy Sersic + broken exponential 

CDd 

14 

2 

Sersic + double exponential 

BDD 

15 

3 

Double Sersic + exponential 

BSD 

16 

3 

Triple Sersic 

BSS 

17 

3 

Double Sersic + broken exponential 

BSDd 

18 

3 


Table 2. Table of the multi-component models used during decomposi¬ 
tion analysis, including the number of independent structural components, 
ricomp, and number of free parameters, k. 


ter galaxies derived from a wide variety of independent imaging 
sources (MegaCam, HST ACS, SDSS, Pan-STARRS, etc.). 


3 ANALYSIS 
3.1 Decomposition 

Galaxy decomposition was carried out using galfit (version 
3.0.4; [Peng et al.||2010| ) with an automated python wrap-around 
derived from agonii (Automated Galfitting of Optical and Near 
Infra-red Imaging; Paper I). Details of this fitting procedure (in¬ 
cluding description of the extraction and calibration of input Mega¬ 
Cam data products) can be found in Appendix |B2|and Paper I. 
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Figure 1. Diversity of surface brightness profiles observed in the Coma cluster galaxies. Six representative galaxies are shown. The first panel of each row is 
the galaxy thumbnail image contour plot (100'' x 100"). The second panel shows the areas masked in the fitting procedure. The ellipses in these two panels 
denote the extent of the galaxy area fitted. The third panel shows the major axis surface brightness profiles (30° width wedges) and errors as derived from 
the stack MegaCam images; labelled by the SDSS Ob j ID number. The measured profiles derived from the individual MegaCam exposures are shown as grey 
lines. The one per cent level of the sky brightness is shown as the horizontal dashed line. The vertical grey bai' shows the radius limit used in the fitting. 
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In the present work, we fit galaxies with a range of analyti¬ 
cal models in order to thoroughly explore the diversity of internal 
galaxy structures. These candidate models are comprised of 1 to 
3 structural components, each described by one of four functional 
forms (see Table^: exponential ‘discs’ (‘D’), general Sersics (‘5’), 
boxy Sersics (‘C’), and broken discs (‘Dd’). Note that the central 
5 component in any model is referred to as the bulge, and labeled 
as ‘6’. Conversely, non-central 5 components in models containing 
a disc are referred to here as ‘bars’. However, this convention does 
not explicitly require a stellar bar structure. As such, a ‘bar’ may 
also correspond to a lens or oval structure. The broken disc com¬ 
ponent (see Appendix |B4| ( comprises inner and outer exponential 
discs (with differing scale lengths, iis.in and i?s,out) connected by 
a smooth transition at a break radius (rurk). 

The ten candidate multi-component models considered in this 
work are catalogued in Table Sersic-only (hereafter ‘5”) and 
bulge -I- disc (hereafter ‘BD’) models are unchanged from those pre¬ 
sented in Paper I. In addition, we present the decomposition results 
when boxy bulge -l- disc (hereafter "CD’), double Sersic (hereafter 
‘6S’), bulge -I- double disc (hereafter "BDD’), bulge -l- bar -l- disc 
(hereafter ‘BSD’), bulge -l- double Sersic (hereafter ‘BSS’) models 
are also considered. Three further models variants implement the 
‘broken disc’ profile: bulge -l- broken disc (hereafter ‘BDd’), boxy 
bulge -I- broken disc (hereafter ‘CDd’) and bulge -r bar -l- broken 
disc (hereafter ‘BSDd’). In order to avoid fitting bias due to the 
choice of initial parameter values, model inputs are based on the 
best fit parameters of simpler model types (e.g. BSD input derived 
from best BD fit). This iterative build up of model complexity sig¬ 
nificantly improves reliability of the measured galaxy properties, 
particularly for highly-degenerate multi-component models. 

We use [Alien et al.|p006^ types to describe the relationship 
between the bulge (i.e. innermost Sersic) and (exponential or bro¬ 
ken) disc profiles (see Appendix |B3^ . This convention is also used 
for 3-component model systems, as fewer constraints are place on 
bar/disc or bulge/bar morphology. The only exception in which at¬ 
tention is paid to these profile interactions is where profile inver¬ 
sion implies incorrect interpretation of the model components (i.e. 
‘Type 4’ bulge/bar or bar/disc structures; see also Appendix |B5[ ). 

3.2 Model Selection and Results Filtering 

Sample filtering is applied to the fitting results (similar to Paper I) 
in order to isolate a sample of accurately-fit galaxies. A key step in 
this process is the selection of best-fit models which are meaning¬ 
ful descriptions of each galaxy’s underlying morphological struc¬ 
ture, ensuring that all structural components are statistically justi¬ 
fied. Galaxy models are assessed on both goodness-of-fit (i.e. en¬ 
suring that a galaxy is neither under-fitted nor over-fitted), and suit¬ 
ability of component structures (i.e. rejecting components with un¬ 
realistic parameters, or which do not measure the intended target 
substructure). By removing such instances of dissonance between 
the galaxy and model stuctures, the reliability of multi-component 
analysis results is vastly improved. A detailed description of the 
galaxy filtering conditions, and a flow chart illustrating the overall 
filtering process is presented in Appendix |B5| 

Galaxies are initially assessed for asymmetry (via the A pa- 
rameter; [Homeier et al.|20()6t and contamination (via image mask 
fraction, /mask) to ensure robust measurements of galaxy prop¬ 
erties. Highly asymmetric galaxies, or galaxies strongly contam¬ 
inated by neighbouring sources cannot be reliably fit by smooth, 
symmetric models, and are thus removed from consideration. Due 
to high parameter uncertainty, galaxies are also removed if their 


best-fit models are poorly-fitted (high xS). highly inclined to the 
line of sight (from the axis ratio of the outer component), or if a 
model component contributes less than 5% of the total galaxy lumi¬ 
nosity (component fraction, C/T< 0.05). Additional filtering con¬ 
ditions are placed on broken disc galaxies to ensure that both the 
inner and outer disc contribute significantly to the overall galaxy 
profile, and to avoid erroneous regions of parameter space. 

Selection between alternative candidate models is made using 
the Bayesian Information Criterion (BIG; |Schwarz||1978^ , calcu¬ 
lated over independent resolution elements (see details in Paper I). 
This is defined as 

where is the standard (un-reduced) fitting chi-squared, k is 
the number of model parameters (degrees of freedom), ripix is the 
number of image pixels used during fitting, and Apsf is the area 
of a resolution element (in pixels). Here, Apaf is calculated as the 
area within two standard deviations (a) of the psf image centre, as 
measured by fitting a Gaussian model. For a set of candidate mod¬ 
els, the model with the lowest BICres maximises goodness-of-fit 
without introducing unnecessary free parameters (hereafter ‘best- 
fit’ model). This ensures that each of the best-fit model components 
provide a statistically significant improvement to x^- 

Measurement error in Apsf leads to an associated uncertainty 
in BICres (cTres). Therefore, a 3crres reduction in BICres is required 
before a more complex (higher k) model is accepted as a statistical 
improvement over a simpler model. This 3a selection condition is 
based on comparison with by-eye classification, and is discussed in 
further detail in Paper I, and |Head] ( |2014| l. Here, (Ties is based on the 
scatter in Apaf (typically ~ 3%), as measured across multiple star 
images. For example, a galaxy fit by S and BD models yields BICres 
values of 3500 and 3450 (respectively) with an associated (Jres of 
10. Since ABICres = 50 > Sures, the addition of the exponential 
disc component is a statistically significant improvement to the fit, 
and hence measures a distinct photometric structure. 

4 RESULTS 
4.1 Best-fit Models 

A wide mix of best-fit model morphologies are found for the 631 
galaxies (570 Coma sample -l- 61 blue Coma galaxies) investigated 
(see example plots in Appendix [A^. The fractions of galaxies best 
described by each candidate multi-component model are illustrated 
in Figure]^ From the initial sample (N — 631), 162 are best fitted 
by a Sersic-only model (26%), 102 are best fitted by BD (16%), 43 
are best fitted by BS (7%), three are best fitted by CD (< 1%), 18 
are best fitted by BDd (3%), three are best fitted by CDd (< 1%), 
43 are best fitted by BDD (7%), 97 are best fitted by BSD (15%), 
98 are best fitted by BSS (16%), and 62 are best fitted by BSDd 
(10%). Thus, the majority of Coma cluster galaxies (58%) have 
morphologies more complex than the simple 5 and BD models. 

Many of these complex structure galaxies were considered to 
be (archetypal SO) bulge -l- disc systems in Paper I. In total, 51 
‘archetypal’ galaxies (from N = 200; 25.5%) remain best-fit by 
a BD model, while 129 (64.5%) require more complex models and 
20 (10.0%) are demoted to a single-Sersic model (due to the more 
stringent model selection tests in the present work). The fractions of 
prior ‘atypical bulge -l- discs’ adequately fit by a BD model (21.1% 


© 2015 RAS, MNRAS 000.[T|j26l 










6 J.T.C.G. Head et al. 



Figure 2. Histogram of best-fit model type fractions for the initial (N = 631) and filtered {N = 478) samples. Error bars are 68% confidence limits. 


of 137), requiring more complex models (73.0%), and demoted to 
S (5.8%) are similar. Note that many galaxies classed as ‘unstable’ 
previously are best-fit here by complex 3-component models (68% 
of 128). This is due to significant reductions in as additional 
structural components are accounted for. 

After sample filtering (see Section [T^ , 478 galaxies from the 
sample of 631 (76%) remain. The 153 galaxies removed by filter¬ 
ing comprise: 80 galaxies removed due to asymmetry or contami¬ 
nation, 23 galaxies with high 13 highly-inclined galaxies, three 
galaxies with anomalous outer discs due to Rs,oMt > 0.1 rurk, two 
galaxies with anomalous inner discs due to rurk < 5", and 32 
galaxies removed due to inverted Sersic/disc components (i.e. disc- 
dominated at low radii, Sersic-dominated at large radii). From the 
remaining filtered sample, 134 galaxies are best fitted by a Sersic- 
only model (28.0%), 94 are best fitted by BD (20%), 34 are best 
fitted by BS (7%), two are best fitted by CD (< 1%), 13 are best 
fitted by BDd (3%), none are best fitted by CDd, 38 are best fitted 
by BDD (8%), 54 are best fitted by BSD (11%), 70 are best fitted 
by BSS (15%), and 39 are best fitted by BSDd (8%). 


Note that in total, 93 galaxies (20 ± 2%) are well-described 
by ‘barred’ models (BSD, BSDd). This barred fraction for Coma 
is sign ificantly lower than the value reported in |Lansbury et al.| 
(2014 1 from either decomposition (72)'(g%) or ellipse (48 ± 6%) 
analyses. This difference cannot be reconciled, even if BSS models 
are included in the ‘barred’ sample (yielding 34^2% bar fraction). 
Flowever, if the present sample is restricted to only contain galaxies 
with D80 morphological classifications (as in |Lansbury et al.|2014| (, 
then the barred fraction (including BSS galaxies) rises to 63 ± 4%. 
This fraction rises further if only D80 SOs (including SO/a, E/SO) 
galaxies are considered, yielding bars in 71 ± 5% of galaxies. As 
the D80 catalogue only covers the bright end of the Coma sample 
(Mg < —18), the bar fraction increase for D80 galaxies indicates a 
significantly decreasing bar detection rate for faint galaxies. How¬ 
ever, the lower bar detection rate relative to |Lansbury et al.| ( |2014| l, 
particularly if BSS galaxies are not considered ‘barred’, reflects the 
more stringent conditions for accepting a more complex model in 
the present work. 


Structural biases in decomposition studies with overly- 
simplistic galaxy models can be quantified by artificially limiting 
the range of candidate models considered during model selection. 


If model selection were repeated without considering double/triple 
Sersic models, 94 out of 631 galaxies (20%) would be identified as 
best fitted by models including a broken disc (BDd, CDd, BSDd). 
Of these, 42 galaxies are better dsecribed by a double or triple 
Sersic model. Thus, a 2D decomposition analysis falsely reports 
broken disc models 45% of the time if only models with exponen¬ 
tial discs are considered. Alternatively, if 3-component models are 
excluded from consideration for 2D analysis, a 2-component model 
is preferred in 342 galaxies (72% of total). However, 199 of these 
galaxies would be better fit by a 3-component modej^ Thus, 2D 
analysis selects an overly-simplistic 2-component model 58% of 
the time if 3-component models are not considered. 

Of all 141 galaxies with 2-component structures, 82 galax¬ 
ies (58%) exhibit Type 1 (i.e. ‘archetypal’ inner -l- outer compo¬ 
nent) profiles, while 44 Type 3 profile galaxies (recurrent bulge; 
31%) make up the second most common structural type. For 3- 
component systems (N = 201), a larger proportion of galaxies 
(130; 65%) are characterised by Type 1 bulge/disc structures, while 
only 55 galaxies (27%) had Type 3 bulge/discs. Thus, 2- and 3- 
component galaxies have archetypal bulge -l- disc structures in the 
majority of cases, with recurrent bulges (dominant over their discs 
at large radii) being the second most common structure. 

If a BD model is forced on the 478 galaxies in the filtered sam¬ 
ple, 214 (62%) yield a Type 1 (archetypal) profile, and 72 galax¬ 
ies (21%) correspond to a Type 3 (recurrent bulge) profile. These 
profile fractions do not change significantly if galaxies best-fit by 
a 2 or 3-component model are considered separately (Type 1/3: 
66%/20% for 2-component galaxies, 59%/22% for 3-component 
galaxies). Thus, Type 3 BD profiles do not intrinsically represent 
underfit galaxy structures, but rather a structural morphology dis¬ 
tinct from archetypal bulge -l- disc systems. 

In summary, thorough 2D decomposition analysis reveals a 
rich range of galaxy structures in the Coma cluster, with 3(-l-) struc¬ 
tural components required in 42% of galaxies. 2-1- component struc¬ 
ture systems were well-represented by archetypal (central) bulge -l- 
(outer) disc morphologies in the majority of cases (N = 202), 


^ Disparity between this value and the 201 3-component models reported 
above is due to two galaxies which would be better fit by a single Sersic 
model, if 3-component models are excluded. 
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Model 

Parameter 

Comp. 1 

Comp. 2 

Comp. 3 

S 

n 

1.90 ±0.05 

- 

- 


Re [kpc] 

1.99 ±0.07 

- 

- 

N = 134 

<? 

0.63 ±0.02 

- 

- 

rui = 16.58 ± 0.03 

C/T 

1.0 

- 

- 

BS 

n 

2.12 ±0.30 

0.66 ± 0.10 

- 


Re [kpc] 

5.27 ± 0.90 

14.57 ± 0.90 

- 

II 

CO 

<] 

0.66 ±0.04 

0.79 ±0.04 

- 

rrii = 15.76 ±0.12 

C/T 

0.35 ±0.04 

0.65 ±0.04 

- 

BSS 

n 

2.17 ± 0.23 

0.43 ±0.13 

0.55 ±0.04 


Re [kpc] 

4.20 ±0.55 

12.68 ±0.76 

25.80 ± 1.75 

II 

0 

<? 

0.71 ± 0.03 

0.48 ± 0.04 

0.58 ±0.04 

rrii = 14.45 ± 0.08 

C/T 

0.42 ±0.02 

0.26 ±0.02 

0.27 ± 0.02 


Table 3. Table of the median structural parameter values for multi-Sersic model galaxies (S, BS, BSS), indicating the half-light radii, Sersic indices, component 
axis ratios (q) and component fraction (C/T) of each model component. In addition, the median total apparent magnitude (m;), and number of galaxies (N) 
are given for each model type. 


including 52 galaxies which exhibited broken disc profiles. This 
broken disc fraction would be overestimated, however, if multi- 
component models (including double/triple Sersic systems) were 
not considered during model selection. 


4.2 Multi-Sersic structures 

From the initial sample of 478 galaxies (filtered to exclude badly-fit 
galaxies; see Section [T^ , ~ 50% were best described by a model 
comprising one or more Sersic components (28% 5; 7% S5; 15% 
B55lj In this section, we briefly discuss the structural results for 
these multi-component Sersic galaxies (Table[^. Note that average 
galaxy luminosity increases with number of model components, 
highlighting the strong S/N and spatial size dependence of multi- 
component structure detection. 

Figurej^illustrates the distributions of n for inner (left panel) 
and outer (right panel) Sersic components. In all Sersic model vari¬ 
ants, the central structure is compact and has a ‘pseudobulge-like’ 
(n ~ 2) profile. For S model galaxies, the Sersic structure is equiv¬ 
alent to a ‘naked’ bulge for BD galaxies, albeit a factor of 2 x larger 
(Re ~ 2 kpc). The central ‘bulges’ of both BS and BSS galaxies are 
consistent in size, but larger on average than a single Sersic (Re ~ 4 
kpc). Note that very few inner structures refer to a classic n = 4 
(de Vaucouleur’s) profile. If BSS (or other 3-component) galaxies 
were force fit by a single Sersic structure, however, the resulting n 
distribution would extend to n ~ 8 , peaking strongly for n = 3-4. 
Hence, de Vaucouleur’s profile may arise from underfitting more 
structurally complex systems. 

Outer Sersic structures have Gaussian-like profiles (n ~0.5) 
on average, although a weak tail exists in the n distribution to¬ 
wards higher values (Figure right panel). Since an outer com¬ 
ponent with n = 1 would be described by a BD model, the 
1.00 < n < 1.25 bin is empty for outer structures. If the disc 
n is allowed to vary for these bulge -H disc galaxies (i.e. fitting a 
BS model), then a continuous distribution of outer structure Sersic 
index becomes apparent (red dashed histogram in Figure]^. The 
resulting ‘disc’ n distribution covers the range 0.5 < n < 1.5, 


but peaks strongly at n = 1 (median value: 1 . 00 , standard devia¬ 
tion: 0.24f] Hence, a subset of the outer structures considered in 
this section may represent (n 7 ^ 1) discs. This is supported by the 
detection of rapid gal axy rotation (14 > 100 km s ^) for 11 of 


our 70 BSS galaxies in 


Rawle et al. 


(2013 


I of their SO sample). 


However as a practical choice, only n = 1 discs will be considered 
in discussion of disc structure in later sections due to the uncertain 
nature of outer Sersic structures. 

The outer and middle Sersic structures of BS and BSS galaxies 
are both > 10 kpc larger than ‘bulges’, but represent drastically dif¬ 
ferent fractions of their parent galaxy’s total luminosity (~ | and 
~ j respectively). Conversely, the outer structure of BSS galaxies is 
comparable in luminosity to the middle Sersic , but is an additional 
10 kpc larger. As such, BSS galaxies are structurally equivalent to 
BS galaxies with the addition of an outer Sersic structure. The outer 
Sersic structures may be the remnants of past merger events. As 
such, the distinction between BS and BSS may be a difference in 
the number of major merger events experienced in the past. 

By comparison, the triple Sersic structures measured by 
|Huang et al.| ( |2013] H13) in a small sample of nearby (visually- 
selected) ellipticals consist of a faint, compact central object (Re < 
1 kpc), a middle component (Re ~ 2.5 kpc), and a dominant outer 
envelope (Re ~ 10 kpc). If the compact components are neglected, 
the HI3 structures are comparable with the multi-Sersic models in 
the present work, albeit with smaller bulges, and more centrally- 
concentrated outer profiles (n ~ 1-2 in H13). This discrepancy 
in outer component n may indicate that the outer profiles in H13 
encompass multiple distinct Gaussian structures (e.g. both outer 
Sersic components in our BSS galaxies). Alternatively, given the 
low local environment density of galaxies in the H13 sample, the 
increased detection rate of weak additional outer Sersic structures 
in the present work may instead reflect a more active merger his¬ 
tory of present-day Coma cluster galaxies. This is supported by the 
higher average bulge size in the present work, as mergers will also 
increase bulge Re- 

If the BS and BSS models considered in this section repre- 


^ Note that outer Sersic n 7 ^ 1 for these models. Thus, we exclude the BD 
and BSD models as special cases of BS and BSS. 


^ Recall however, that these changes in outer profile n do not yield sta¬ 
tistically significant improvements to the goodness-of-fit relative to fixing 
n = 1 . 
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Figure 3. Histograms of Sersic index, n, for multi-component Sersic galaxies. Left: the n distribution for inner-dominant Sersic structures (S, BS, BSS), 
divided by best-fit model type. Right: the n distribution of outer-dominant Sersic structures (BS, BSS) divided by best-fit model type. The equivalent distribu¬ 
tions including middle Sersic structures (BSS) or outer Sersic n for BS fits to galaxies best-fit by a BD model are included as dashed blue and red histograms 
respectively. 


sent the multi-component structures of traditional elliptical galax¬ 
ies, then such galaxies comprise a (relatively) compact pseudob¬ 
ulge (~ 5 kpc) around which large (~ 10-20 kpc) outer (Gaussian) 
structures have been assembled. These compact central structures 
are 2-3 x larger than ‘red nugget’ objects (~ 1-2 kpc) detected at 
high redshift (z ~ 2; |DamJanov et al.|[2009| l. Thus, if the multi- 
component Sersic galaxies observed in Coma in the present work 
evolved from red nuggets, then their bulge structures must have 
experienced significant size growth (‘puffing up’). However, the 
total effective radii for BS and BSS galaxies (estimated from the 
combined luminosities of all model components, assuming align¬ 
ment of component PAs) is ~ 10-11 kpc on averag^ suggesting 
an even more drastic growth mechanism (~ 6 x, consistent with 
[van Dokkum et al.|2014t . 

In summary, galaxies comprising multiple Sersic structures 
(with outer n 7 ^ 1 ) resemble a compact central psuedobulge (rem¬ 
iniscent of single Sersic systems; n ~ 2, ~ 4 kpc) embedded 

in extended Gaussian (n ~ 0.5) envelopes. The combined effective 
half-light radii of these systems typically exceeds 10 kpc. Thus, if 
multi-Sersic galaxies evolved from compact ‘red nuggets’ as de¬ 
tected at high redshift, then such systems must have experienced a 
> 6 X increase in size. 


4.3 Freeman Disc Type Fractions 

Galaxies with single disc-like outer profiles (BD, CD, BDd, CDd, 
BSD, BSDd) were categorised by their disc types (i.e. Freeman 
Type I, II, or III). In total, 202 valid disc galaxies are contained 
within the sample after filtering. Of these, 150 galaxies (74 ± 3%) 
have Type I (untruncated) discs, 25 galaxies ( 12 l 2 %) have Type 
II (truncated) discs, and 27 galaxies have Type III (anti¬ 

truncated) discs (Figure]^ left panel). Compared to the disc type 
fractions reported in the Virgo cluster (Type I: 46 ± 10%, Type 


® Even if the outermost structures in BSS galaxies were dismissed as fitting 
artefacts, the total i?e of such systems would remain in excess of 7 kpc. 


II: 01q%, Type IIF 54 ± 10%;|Erwin et al.|2012i, we detect sig¬ 
nificantly more Type I and II discs in Coma, but fewer Type III 
discs. By comparison, the field SO sample in |Erwin et al.|pbl2| > 
yields significantly fewer Type I discs (26)'lg%), but greater Type 
II ( 2810 %) and III (46 ± 7%) fractions than the Coma sample. 


If considered separately. Type I discs are found more fre¬ 
quently in unbarred (BD, CD, BDd, CDd', 88 ± 3%) galaxies than 
those containing bars (BSD, BSDd', 58 ±5%). Consequently, barred 
galaxies have a greater fraction of Type II and III discs ( 23^^% 
and 19 ± 4%) than galaxies without bars (4l^% and 8 l 2 %)-|Erwin 
|et al.| also reported a decrease in the Type I fraction for barred Virgo 
galaxies (23^7%), however the increased barred Type II fraction 
in this work only widens the disparity between Coma and Virgo 
Type II disc detection. Note that no strong correlation (Pearson’s 
p ~ 0.3) is detected between bar and broken disc axis ratios (q) of 
BSDd galaxies, indicating that these model components are struc¬ 
turally distinct. Thus, the detection of a large number of broken 
discs in barred galaxies is not an artefact of of overfitting (i.e. via 
coupling of the inner disc to the Sersic bar profile). 


To test the variation of disc type with environment, the fil¬ 
tered Coma sample was sub-divided into core, intermediate, and 
outskirts samples based on galaxy distance from the cluster centre 
(fcluster ^ 0.8 MpC, 0.8 ^ nduster ^ 1.6 MpC, and ^cluster ^ 1.6 
Mpc respectively; Figure]^ right panel). These clustercentric radial 
ranges are selected such that each sample has approximately equal 
occupancy (N = 70, 68, and 64). In all three samples. Type I discs 
form the vast majority, with a slightly increased Type I disc frac¬ 
tion for outskirt galaxies (73lg%, 72j')g%, and 78lg% for core, 
intermediate, and outskirt galaxies). Type II and III disc fractions 
are consistent across all radial samples (Type II: 10 ^ 3 %, 

Type III: 17l^%, 12t*%, although slight peaks 

in Type II and Type III disc fractions are apparent in the intermedi¬ 
ate and core samples (respectively). 


In summary, greater fractions of Freeman Type I (untruncated; 
74%) and Type II (truncated; 12%) discs were detected in the 
present work than have been reported previously in the Virgo clus¬ 
ter. Conversely, the measured fraction of Type III (anti-truncated) 
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Type II Type III 



r [arcsec] r [arcsec] 


Figure 4. Example (major axis) surface brightness profiles for six galaxies (3 BDd, 3 BSDd) best fitted by broken disc models (3 Type II, 3 Type III). 
Upper panels: The i band surface brightness as measured from the galaxy thumbnail (black points) in wedges of elliptical annuli. The corresponding model 
components are indicated as solid lines (black/red/green/blue/cyan: total/bulge/bar/inner disc/outer disc). Lower panels: Model residual (in counts). Error 
bars in both plots (gray) are the standai'd error on the mean surface brightness in each wedge. All examples include the ratio of the outer and inner disc scale 
lengths (-Rs,out/Hs,in). 
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Figure 5. Histogram of outer disc type fractions for Left: banned (BSD, BSDd) and unbarred (BD, BDd) bulge + disc models, Right: galaxies in the core 
(fcluster < 800 kpc), intermediate (0.8 < fciuster < 1-6 Mpc), and outskirts (rduster > 1-6 Mpc) cluster samples. Types I, II, and III refer to unbroken, 
truncated, and anti-truncated discs respectively. Error bars are 68% confidence limits. 


discs in Coma (13%) was lower than Virgo. The majority of galax¬ 
ies with Type II or III discs also contain galaxy bars (Type II: 89%; 
Type III: 71%), compared to less than half of galaxies with unbro¬ 
ken discs (Type I: 42%). No significant variation in Type I/II/III 
fraction was detected with local environment within Coma. 

4.4 Freeman Type I, II, and III galaxy structures 

In this section, the galaxy sample is divided by Freeman type to 
investigate differences in internal structure for galaxies with un¬ 
truncated (Type I), truncated (Type II), or anti-truncated (Type III) 
discs. We consider all models with a single (exponential) disc com¬ 
ponent which dominates (relative to the bulge) at large galaxy¬ 
centric radii (BD, CD, BDd, BSD, BSDd). Here, the distributions 
and trends in structural parameters for galaxies with (single) disc- 
dominated outer regions (BD, CD, BDd, BSD, BSDd) are investi¬ 
gated. In order to ensure that model parameters are measured from 
consistent structures (i.e. exponential/broken exponential compo¬ 
nents measure galaxy disc properties), we only consider galaxies 
with archetypal bulge/disc models (i.e. Type 1; central bulge -l- 
outer disc). This reduces the sample of analysed galaxies to 146 
(67 2-component galaxies, 79 3-component galaxies), of which 97 
galaxies have Type I discs, 18 have Type II discs, and 24 have Type 
III discs. The average structural properties and total magnitudes of 
archetypal galaxies containing discs of each type are summarised in 
Table Example surface brightness profiles for (BDd and BSDd) 
galaxies with Type II and Type III discs are presented in Figure 
1^ As a convenient shorthand, we hereafter use the phrase ‘Type 
FII/III galaxy’ to refer to galaxies containing Freeman Type I/II/III 
discs. 

4.4.1 Central components of Type I/II/III galaxies 

The bulge and bar Sersic indices for galaxies of each disc type are 
presented in Figure]^ Bulge n is smaller (on average) in Type I 
(1.89±0.08) galaxies those with Type II (2.32±0.14) broken discs, 
and consistent with the bulges of Type III galaxies (1.74 ± 0.19). 
By comparison, bar Sersic index is consistent across all Freeman 
types (0.44 ± 0.04 for Type I, 0.44 ± 0.04 for Type II, and 0.43 ± 


0.03 for Type III). Thus, while consistent bar profiles are measured 
independent of disc type, the bulge profile shape depends on disc 
structure. Note that the Type I averages are calculated from galaxies 
in the magnitude range—19 < Mi < —22 for consistency with the 
range of Type II and III galaxy luminosities. 

With increasing galaxy luminosity, no significant variation in 
bar n is detected for any galaxy type. However, the bulges of both 
Type I and Type III galaxies become more centrally-concentrated 
(higher n) for more luminous galaxies. Similar n-Mi slopes are 
measured for both galaxy types (consistent with the equivalent 
trend measured previously for archetypal BD models in Paper 
I). The reverse trend (lower n for higher galaxy luminosity) is 
measured for Type II galaxies. While this trend is not significant 
(~ 1.5a), it remains discrepant with the measured Type I/III trends 
at a 3a level. Thus, the bulges of galaxies with truncated discs 
are structurally distinct from those found in galaxies with untrun¬ 
cated, or anti-truncated discs. This is analogous to the distinct n- 
luminosity trends measured in the previous section for barred and 
unbarred galaxies. However, as barred galaxies comprise approxi¬ 
mately equal numbers of Type II (231®%) and III discs (19 ± 4%), 
this apparent bulge n bimodality is not strongly related to the pres¬ 
ence of a bar component. 

Half light radii for the bulges and bars of Type I, II, and III 
galaxies are presented in Figure The bulges of Type I and III 
galaxies show no significant size difference on average (Type I: 
0.57 ± 0.05 kpc. Type III: 0.54 ± 0.05 kpc), while Type II galaxies 
have systematically larger bulges (0.96 ± 0.07 kpc). No notable 
trends in bulge size with galaxy luminosity is noted for galaxies of 
any Freeman type. 

The bars in Type II galaxies are systematically larger on 
average (2.91 ± 0.22 kpc) than those found in Type I galaxies 
(1.95 ± 0.18 kpc), but similar to the bars of Type III galaxies 
(2.82 ± 0.49 kpc). As with bulge components, no significant size- 
luminosity trends are noted for galaxy bars. Thus, large galaxy bars 
are found more frequently in galaxies with broken (truncated/anti¬ 
truncated) discs, regardless of total galaxy luminosity. 

In summary, the bulges of galaxies with Type II discs have sys¬ 
tematically larger n and Re than the bulges of Type I or III galaxies. 
In addition, no significant bulge n-luminosity trend is detected for 
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Disc Type 

Parameter 

B 

S 

D/Dd 

Type I 

n 

1.52 ±0.09 

0.44 ± 0.04 

1.0 


Re [kpc] 

0.62 ±0.04 

1.94 ± 0.16 

3.68 ±0.18 

N = 97 

<? 

0.65 ±0.02 

0.54 ± 0.04 

0.67 ± 0.02 

mi = 15.82 ±0.10 

C/T 

0.27 ± 0.02 

0.11 ± 0.02 

0.62 ±0.02 

Type II 

n 

2.32 ±0.14 

0.44 ± 0.04 

1.0 


Re [kpc] 

0.96 ±0.07 

2.91 ± 0.22 

9.28 ±1.44 / 4.12 ±0.20 

N = 18 

<? 

0.69 ±0.05 

0.49 ± 0.08 

0.56 ±0.05 

mi = 14.65 ±0.12 

C/T 

0.33 ±0.04 

0.21 ± 0.03 

0.46 ± 0.04 

Type III 

n 

1.74 ± 0.19 

0.43 ± 0.03 

1.0 


Re [kpc] 

0.54 ± 0.05 

2.82 ± 0.49 

3.58 ±0.20 / 5.62 ±0.34 

N = 24 


0.70 ±0.04 

0.39 ± 0.07 

0.54 ± 0.05 

mi = 14.49 ±0.17 

err 

0.24 ± 0.04 

0.14 ± 0.02 

0.62 ±0.03 


Table 4. Table of the average structural parameter values for Type I/II/III archetypal disc galaxies with simple exponential discs {BD, CD, BSD) and broken 
exponential discs {BDd, BSDd). The average Sersic indices (n), half-light radii (Re), component axis ratio (q), and component fractions (C/T) are indicated for 
each model component. For Type II and III galaxies, both the inner and outer Re (= 1.678Rs) are included. In addition, the median total apparent magnitude 
(mi), and number of galaxies (N) are given for each disc type. 


Type II galaxies. Thus, the bulges of galaxies with truncated discs 
are distinct in structure and origin from the equivalent components 
in galaxies with either untruncated or anti-tmncated discs. Galaxy 
bars are consistent in profile shape across all Freeman types, but 
have systematically larger Re in galaxies with Type II or III broken 
discs. 

4.4.2 Type I/II/III Disc Properties 

The effective half-light radii for the inner and outer discs (i.e. the 
discs internal and external to the break radius, rurk) of Type II 
(truncated) and III (anti-truncated) galaxies are presented in Fig¬ 
ure with the disc Re for Type I galaxies included in both pan¬ 
els. Note that by definition, the inner Re of Type II/III galaxies 
is larger/smaller than the outer Re, yielding a shallower disc sur¬ 
face brightness profile within/beyond rurk- On average. Type III 
inner discs are consistent in size (3.58 ± 0.20 kpc) with Type I 
discs (3.68 ± 0.18 kpc), and have a consistent size-luminosity re¬ 
lation (despite a ~2x difference in slope). By contrast. Type II 
inner discs are substantially larger (than Type I discs) on average 
(9.28 ± 1.44 kpc), with an extremely steep trend (4.76 ± 3.09 kpc 
per mag) of decreasing inner disc size with increasing galaxy lumi¬ 
nosity, albeit at low significance (~ 1.5cr). 

The outer discs of Type II galaxies have scale lengths (4.12 ± 
0.20 kpc) consistent with Type I discs (for galaxies in the range 
—19 < Mi < —22) on average, while Type III outer discs are 
systematically larger (5.62 ± 0.34 kpc). Outer disc size-luminosity 
relations are similar for both Type II and III discs, yielding size 
increases for more luminous galaxies a factor of approximately 
two times greater than the measured trend for Type I discs. How¬ 
ever, this difference relative to Type I discs is only significant 
(at a ~ 2.5 (t level) for Type II galaxies. The detection of consis¬ 
tent scale lengths (and similar size-luminosity relations) for Type 
I discs. Type II outer discs, and Type III inner discs (in agreement 
with |Laine et al.|2014| l suggests that the outer/inner structures of 
Type II/III discs preserve the structural properties of their progeni¬ 
tor discs. 

The break radius, rurk, is plotted in Figure|^for Type II and III 
discs as a fraction of both inner and outer disc Re. For Type II discs, 
the break radius is a small fraction of the inner disc size (~ 0.25). 


However, since Re is large for these structures, the contribution of 
the (flat) inner disc to the total disc luminosity is non-negligible. 
By comparison. Type III disc break radii are significantly beyond 
the inner disc half-light radius (rbrk ~ 2i?e), indicating that only 
the outer wings of Type III inner disc structures are modified by the 
profile break. Alternatively, both Type II and III profile breaks are 
comparable in size to the outer disc Re. Thus, in both cases, the 
outer structure of broken discs contribute ~ 50% of the light of an 
equivalently-sized untruncated disc. 

In comparison to either disc. Type III breaks occur at smaller 
fractions of disc Re for increasingly luminous galaxies. Note, how¬ 
ever, that this correlation is significant at a > 3a level for outer 
disc Re, but only significant at a ~2.5 ct level for inner disc Re. 
A decreased fractional break radius indicates that a Type III disc 
contains a smaller proportion of the primordial disc. Conversely, 
Type II break radii exhibit a non-significant increase (as fractions of 
both Re) with galaxy luminosity. Thus, the break radius of a Type 
II disc is approximately the same fraction of the inner/outer disc 
size for any galaxy. Note that if Type II discs where Re,in reaches 
the GALFIT limit are excluded, the trend in rurk relative to Re,out 
is made considerably shallower. Hence, rurk increases in size at a 
similar rate to outer disc Re with total galaxy luminosity. 

In summary, the inner discs of anti-truncated (Type III) galax¬ 
ies are consistent in size with the discs of unbroken (Type I) galax¬ 
ies. Conversely, Type III outer discs are systematically larger than 
unbroken discs. Both inner and outer Type III discs exhibit a size- 
luminosity relation consistent with Type I discs. Thus, the inner 
discs of Type III galaxies preserve the properties of the unbroken 
progenitor disc. The inner discs of truncated (Type II) galaxies are 
not consistent in size or size-luminosity trend with unbroken discs. 
This rules out a formation scenario in which physical tmncation 
preserves the primordial disc within the break radius. Conversely, 
the outer discs of Type II galaxies have sizes (and size-luminosity 
relations) consistent with untruncated disc structures. 

4.4.3 Component fractions (C/T) of Type I/II/III galaxies 

In this section, we discuss the component flux fractions (B/T, S/T, 
and D/T for bulges, bars, and discs respectively) for galaxies with 
Type I, II, and III discs. 


© 2015 RAS, MNRAS 000.[T|j26l 












12 


J.T.C.G. Head et al. 



N 

^B/T 

^S/T 

^D/T 

Type I 

56 

-0.01 

± 

0.03 




0.01 

± 

0.03 

2-comp. Type II 

4 

-0.10 

± 

3.65 




0.10 

± 

3.65 

Type III 

7 

-0.12 

± 

0.11 




0.12 

± 

0.11 

Type I 

41 

-0.04 

± 

0.03 

-0.10 

± 

0.03 

0.12 

± 

0.04 

3-comp. Type II 

21 

0.11 

± 

0.07 

0.01 

± 

0.05 

-0.22 

± 

0.06 

Type III 

17 

-0.14 

± 

0.02 

0.07 

± 

0.03 

0.05 

± 

0.03 


Table 5. Table of best-fit component light fraction-luminosity trends C/T per magnitude galaxy luminosity) measured for 2-component (top) and 

3-component (bottom) galaxies. Here, a negative value indicates increasing C/T with luminosity. 


s 

01 

_ai 

OQ 



My [mag] 



-18 


-19 


-20 -21 

My [mag] 


-22 


-23 



My [mag] 



My [mag] 


Figure 6 . Bulge and bar Sersic indices (n) for BD, CD, BDd, CDd, BSD, 
and BSDd model galaxies with Type I (black), II (red), and III (blue) discs 
as a function of total absolute i band model magnitude. Upper plot: Bulge 
Sersic index. Lower plot: Bar Sersic index. Unfilled data points indicate 
flagged galaxies. Large square points are median parameter values in bins 
of Mi, to which a linear trend has been fit. Type I galaxies are indicated by 
small grey points for clarity. 


Figure 7. Bulge and bar effective half-light radii (Re) for BD, CD, BDd, 
CDd, BSD, and BSDd model galaxies with Type I (black), II (red), and III 
(blue) discs as a function of total absolute i band model magnitude. Upper 
plot: Bulge Re. Lower plot: Bar Re. Unfilled data points indicate flagged 
galaxies. Large square points are median parameter values in bins of Mi, 
to which a linear trend has been fit. Type I galaxies are indicated by small 
grey points for claiity. 


Measured across all (2- and 3-component) galaxies, Type I 
galaxies are strongly disc-dominated (median D/T = 0.62 zb 0.02), 
with (sub-dominant) bulges (B/T = 0.28 zb 0.02) and weak bar 
components (S/T = 0.10 zb 0.01). The corresponding component 
fractions for Type III galaxies are measured to be consistent with 


Type I galaxies (D/T= 0.62 ± 0.03, B/T = 0.24 zb 0.02, S/T = 
0.14 zb 0.03). By contrast, Type II galaxies have a diminished disc 
light fraction on average (D/T= 0.46 zb 0.04), with corresponding 
increases in bulge (B/T= 0.33 zb 0.04) and bar (S/T= 0.21 zb 0.03) 
fractions. Note, however, that these averages are heavily biased by 
the lack of a bar (i.e. S/T= 0.00) in 2-component galaxies. If the 
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Figure 8. Inner and outer disc effective half-light radii (Re) for BD, CD, 
BDd, CDd, BSD, and BSDd model galaxies with Type I (black), II (red), 
and III (blue) discs as a function of total absolute i band model magnitude. 
Upper plot: Inner disc R,, . Lower plot: Outer disc R^ . Unfilled data points 
indicate flagged galaxies. Large square points are median parameter values 
in bins of Mi, to which a linear trend has been fit. Type I galaxies ai'e 
indicated by small grey points for clarity. 


average is calculated from only 3-component galaxies, then bar 
light fraction (S/T) increases significantly for all three disc types 
(Type I: 0.24 ± 0.02; Type II: 0.25 ± 0.02; Type III: 0.20 ± 0.03). 
The corresponding disc light fractions (D/T) decrease on average 
for 3-component galaxies (Type I: 0.51±0.03; Type II: 0.43±0.04; 
Type III: 0.55 ± 0.03), while average bulge fractions (B/T) are not 
significantly changed (Type I: 0.25 ± 0.02; Type II: 0.31 ± 0.04; 
Type III: 0.20 ± 0.03). 

The best-fit component light fraction trends with galaxy lumi¬ 
nosity (2-component and 3-component galaxies considered sepa¬ 
rately) are presented in Tablej^ For 2-component galaxies, no sig¬ 
nificant trends are noted in Type I galaxy B/T or D/T, while no 
conclusions can be drawn for Type II and III galaxies due to small 
sample sizes. However, with increasing luminosity, 3-component 
Type I galaxies become significantly more bar-dominated (—0.10± 
0.03), and less disc-dominated (0.12 ± 0.04). Conversely, 3- 
component Type II galaxy disc light fraction and Type III galaxy 
bulge light fraction increases with galaxy luminosity (—0.22±0.06 
and —0.14 ± 0.02 respectively). 



M, [mag] 

Figure 9. Break radius, as a function of the total absolute i band 

model magnitude for galaxies with Type II (red) and III (blue) disc models 
relative to the inner (top) and outer (bottom) disc iJe. Unfilled data points 
indicate flagged galaxies. Large square points are median parameter values 
in bins of Mi, to which a linear trend has been fit. 


These component light fraction-luminosity trends can be used 
to estimate whether the distinction between faint and bright galax¬ 
ies is dominated by the luminosity difference of one particular com¬ 
ponent. This can characterise, for example, whether the difference 
between an average galaxy and an equivalent galaxy one magnitude 
brighter is primarily due to an increase in bulge or disc luminosity. 
Hence, we will determine whether the apparent differences in C/T 
trends between Freeman Types corresponds to intrinsically differ¬ 
ent component light scaling relations. 

For two galaxies separated in total galaxy luminosity by one 
magnitude (Mo — M — 1.0), the fractional difference in the lu¬ 
minosity of a particular component (C) can be parametrised as 
xc = Lc ILc,q- For example, if the galaxy luminosity difference 
in Type I galaxies at Mi,o = —20 and Mi = —21 was caused by 
the bulge and disc components being 3x brighter at Mi = —21 
(but bars being as luminous in both cases), then xb = ‘A, xs = 1, 
and xd = 3. 

The reported C/T slopes (Ac/ t; Table|^ can be expressed as: 
Ac/t = Co/To - C/T = Co/Tq (l " fj) - (2) 
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N 

xbIxt 

xsIxt 

xdIxt 

Type I 

41 

1.16 it 0.12 

1.41 it 0.13 

0.76 it 0.08 

Type II 

21 

0.65 it 0.23 

0.96 it 0.20 

1.51 it 0.15 

Type III 

17 

1.56 it 0.09 

0.66 it 0.16 

0.91 it 0.06 


Table 6. Table of approximate fractional component luminosity changes for 
3-component Type I, II, and III galaxies as total galaxy luminosity increases 
(XT = 2.5). 

where fractional difference in total luminosity, xt = 2.5 across 
one magnitude. Note that xc jxT i^s, greater than unity if component 
luminosity increases at a greater rate than galaxy luminosity. 

Table presents xc values relative to a galaxy of average 
luminosity and C/T (i.e. Cq/To = (C/T)). For a galaxy brighter 
than the average by an arbitrary magnitude difference (Mtot = 
Mtotfi + AMtot), the proportion of the total luminosity difference 
(ALtot; where Ltot.o = Got + ALtot) attributed to each photo¬ 
metric component (ALq) can be estimated using: 

ALc ^ ~ 1) ~ xtAc/t ^2) 

ALtot XT — 1 

The resulting component fractions of the additional galaxy lumi¬ 
nosity is illustrated in Figure[T^for Type I, II and III galaxies. 

For Type I galaxies, bars become more luminous at a sig¬ 
nificantly greater rate than the overall galaxy luminosity (bar lu¬ 
minosity doubles for a 42% increase in total luminosity), while 
discs increase in luminosity at a slower rate (52% increase in disc 
luminosity as galaxy luminosity doubles). However, the luminos¬ 
ity difference for Type I galaxies arbitrarily brighter than the av¬ 
erage is distributed equally between all three structural compo¬ 
nents (from Equation]^ Figure [Toll. For example, an average Type 
I galaxy ({Mi) = —20.3) has a total luminosity of 8.6 x 10® 
Lq. Relative to this average, a 9.6 x 10® Lq Type I galaxy 
(i.e. 10 ® Lq brighter) would have a bulge more luminous by 
(3.3 ± 0.4) X 10® 1 / 0 , a bar more luminous by (3.5 ± 0.3) x 10® 
Lq, and a disc (3.2 ± 0.4) x lO®I /0 more luminous. 

For Type II galaxies, the disc component is the dominant con¬ 
tribution to luminosity growth (70 ± 12 % of ALtot), doubling in 
luminosity for each 32% increase in galaxy luminosity. Hence, the 
disc-total luminosity trend is significantly steeper for Type II discs 
than Type I, indicating a larger difference in disc luminosity be¬ 
tween faint and bright Type II galaxies than for Type I galaxies. 
This implies that fainter Type II galaxies have experienced a greater 
truncation (of light) than intrinsically more luminous galaxies. 

For Type III galaxies, the bulge is the dominant component 
(doubling in luminosity for a 28% increase in global luminosity). 
For an average Type III galaxy, the bulge component’s contribution 
to galaxy luminosity is approximately equal to the (intrinsically 
more luminous) disc. The corresponding bar light contribution is 
minimal (3.7 ± 0.4% of ALtot), indicating approximately equally- 
luminous bars in all Type III galaxies, independent of total galaxy 
luminosity. 

No other components (bulges in Type I galaxies, bulges and 
bars in Type II galaxies, and bars and discs in Type III galaxies) dif¬ 
fer significantly from increasing in luminosity proportional to the 
galaxy (xc ~ 1). Note that since xt = 2.5, no component in Type 
I, II, or III galaxies decreases in luminosity in brighter galaxies. The 
proportions of added galaxy luminosity contributed by each model 
component are illustrated in . 



0.2 


Type I Type II Type Ml 

Figure 10. ‘Additional’ component light fractions (i.e. the proportions of 
the luminosity difference for each component structure per unit galaxy lu¬ 
minosity) for average galaxies with Freeman Type I, II and III disc com¬ 
ponents. Indicates the fraction of light added to the bulge (red), disc (blue) 
or bar (green) per unit galaxy luminosity increase. Type I galaxies - bulge: 
33 it 3%, disc: 35 it 3%, bar: 32 it 4%; Type II galaxies - bulge: 10 it 2%, 
bar: 19 it 2%, disc: 70 it 12%; Type III galaxies - bulge: 50 it 4%, bar: 
4 it 1%, disc: 47 it 6%. 

In summary, 3-component archetypal (central bulge -l- outer- 
dominant disc -I- any bar) galaxies are disc-dominated on average, 
with approximately equal bulge and bar light fractions, indepen¬ 
dent of Freeman disc/galaxy type. The measured trends in compo¬ 
nent light fraction with total magnitude were used to quantify the 
contributions of each structural component to galaxy luminosity. 
All three structural components contribute equally on average to 
the increasing total luminosity in galaxies unbroken discs (Type I). 
However, the bar component exhibits the largest fractional increase 
in luminosity, discs were found to dominate truncated (Type II) 
galaxy luminosities. The corresponding disc-total luminosity trend 
is steeper than for Type I galaxies, which may indicate disc (lumi¬ 
nosity) truncation. Increasing anti-truncated disc (Type III) galaxy 
luminosities correlate strongly with both their bulges and discs. 
Hence, bar luminosity in Type III galaxies is independent of galaxy 
luminosity. 

4.5 Structural Trends with Enviroiunent 

In this section, we investigate variation in the multi-component 
structures of galaxies as a function of the (projected) distance from 
the Coma cluster centre, rduster- Observed Tciuster correlates with 
the time at which a galaxy first entered the cluster environment 
( |Gao et al.|2004[ [Smith et al.|2012]|De Lucia et al.|2012[|Taranu| 

|et al.|20f4l >, albeit with substantial scatter. A cluster-centric radial 
analysis therefore highlights the environment-mediated processes 
that have acted on these multi-component systems, and hence the 
cluster environment’s role in their formation. 

The morphological mix of galaxies varies with position in the 
cluster (Figure [TT|. Neither the fraction of multi-Sersic models, nor 
the fraction of galaxies with broken discs vary significantly with 
cluster-centric radius. Note that, this would not change if Type II 
and III galaxies were considered separately (see Section [43) l. The 
fraction of pure Sersic systems increases towards the cluster centre. 
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Figure 11. The model fractions, /mod (i-6- number of each model type 
total number of galaxies per bin) as a function of radial position in the 
Coma cluster, foiuster ■ Single Sersic objects (S) are plotted in red, galaxies 
with outer exponential discs (BD, CD, BSD) are plotted in blue, galaxies 
with broken discs (BDd, BSDd) are plotted in black, and multi-Sersic sys- 
tems(B5, BSS) are plotted in green. 

while the fraction of (exponential) outer disc galaxies decreases. 
However, neither of these radial morphology trends are significant. 

With increasing distance from the cluster centre, no significant 
variation (> 3a) is detected in the structural properties of Type I, II, 
or II galaxies. However, weakly significant radial trends (~ 2a) are 
detected in (barred) Type I galaxy bar size {Re, 1-3 ± 0.6 kpc per 
f 20 o)> Type II outer disc size (i?e,out; —1.5 ± 0.8 kpc per r 2 oo Q 
and Type III galaxy bulge profile (n; 0.9 ± 0.5 per r 2 oo)- Thus, the 
strucural properties of galaxies with disc-dominated outer regions 
are independent of their radial position in the cluster. Therefore, 
the formation of disc breaks may result from secular, rather than 
environment-mediated processes. 


5 DISCUSSION: DISC BREAK FORMATION 
SCENARIOS 

Here, we briefly investigate the evolutionary origins of archetypal 
broken disc (truncated/antitruncated; Freeman Type II/III) galaxies 
through comparison of their structural and component photometric 
properties to unbroken (Type I) galaxies. As a working hypothe¬ 
sis, we assume a break formation scenario in which Type II and 
III galaxies had Type I discs at some point in the past. While the 
observed present-day Type I discs are not necessarily the progen¬ 
itors of present-day broken discs, all three Freeman type galaxies 
are assumed to have evolved from a common population of pri¬ 
mordial galaxies with (Type I) discs. Thus, characteristics of the 
structural/photometric distributions unique to galaxies of a particu- 


® Note that the sign of this trend indicates increasing outer disc size to¬ 
wards the cluster centre, and hence does not correspond to environmental 
truncation. 
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lar Freeman type can be used to constrain their evolutionary path¬ 
ways. 

The absence of any strong cluster-radial trends in galaxy 
structure disfavours a (cluster) environment-driven origin for disc 
breaks. Furthermore, bar structures appear to be strongly related 
to the formation of Type II and III discs: while one half of all 
Type I galaxies contain a bar (42^5% barred, 58^5% unbarred), the 
bar fraction is considerably higher for Type II (89lg%) and Type 
III (7lt.io%) galaxies. Galaxy bars are also significantly larger 
if their host galaxy has a truncated/anti-truncated disc than if the 
galaxy disc remains unbroken. This implies that either the forma¬ 
tion mechanism induces bar growth, or that bars stabilise discs 
during truncation/anti-truncation, such that the detection of a disc 
break for bright galaxies is more likely if a bar is present. 

The detection in galaxy simulations of inner (and outer) disc 
evolution with time ( [Debattista et al.|2006[ [Minchev et al.|2012[ l 
supports a scenario of stellar (or gas) redistribution. In particu¬ 
lar, the radial angular momentum transfer mechanism proposed in 
|Minchev et al.|pbl2^ would explain the apparent importance of a 
bar component, as such a structure would induce significant gravi¬ 
tational torques in disc gas. The significant increase in bar size for 
more luminous Type II galaxies may therefore suggest a period of 
enhanced star formation in the bar due to gas inflows, or the migra¬ 
tion of disc stars into the bar. 

Systematically larger Type II inner disc scale lengths (and in¬ 
consistent size-luminosity relations) compared to untruncated Type 
I galaxies indicates that the inner discs of Type II galaxies are not 
structures equivalent to Type I discs. This disfavours a scenario in 
which Type II discs represent a truncated system in which the outer 
disc is suppressed relative to the surviving primordial inner disc. 
Furthermore, while D/T is systematically lower in Type II galax¬ 
ies compared to Type I, the fractional change in disc light does not 
dilfer significantly from unity {xo/xt = 0.9 ±0.1 as in Equa¬ 
tion]^ where xt is the galaxy luminosity change between Type I 
and Type II galaxies). Hence, assuming an evolutionary scenario 
in which Type II galaxies evolve from Type I, disc luminosity in¬ 
creases proportional to the ~ 40% increase in total galaxy luminos¬ 
ity (i.e. Type I: mr = 14.8 ±0.1 vs. Type II: mr = 14.5 ± 0.1). 
Intrinsically brighter Type II discs rule out a formation mechanism 
in which Type I discs are physical truncated. This conclusion is 
not compromised by the comparison of present-day truncated and 
untruncated discs unless evolution from primordial to present-day 
Type I galaxies also involves reduction of disc luminosity while 
preserving their untruncated profiles. 

Beyond rurk, Type II discs represent structures reminiscent 
of their primordial Type I discs (see also [Foyle et al.|2008^ . Con¬ 
versely, in the inner region (r < rbrk) disc light has been redis¬ 
tributed such that the profile is flattened relative to a Type I profile. 
Bulges and bars in Type II galaxies are systematically larger than 
those in untruncated galaxies, implying that secular bulge/bar en¬ 
hancement effects are significant for the formation of Type II galax¬ 
ies. Thus, disc stars within rbrk may have been redistributed to form 
a bar and/or grow the galaxy bulge (see jValenzuela & Klypin|2003| ). 
Alternatively, the break formation mechanism may be enhanced via 
interaction with an existing bar, resulting in steeper inner/outer disc 
size trends with luminosity due to the strong bar size - galaxy lu¬ 
minosity relation. 

Consistency in component scale-lengths (and size-luminosity 
trends) between Type III and Type I galaxies implies that the Type 
III inner discs may correspond to undisturbed primordial (Type I) 
discs. Conversely, the (significantly larger) outer disc may represent 
an additional extended structure. Nevertheless, this outer structure 
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maintains a disc-like size-luminosity relation. An evolutionary sce¬ 
nario from Type I (or Type I progenitors) to Type III is supported 
by the consistent bulge and disc component light fractions for both 
disc types, despite Type III galaxies being a factor of 1.9 x brighter 
on average. Hence, bulge and disc luminosities increase propor¬ 
tional to the galaxy luminosity difference between Type I and III 
galaxies. 

Bulge and bar sizes in anti-truncated galaxies are significantly 
larger than those in Type I galaxies, while bulge luminosity in¬ 
creases strongly in more luminous Type III galaxies (xb ~ 1.5xt', 
see Table [^. Thus, similar to Type IIs, the formation of Type III 
galaxies involves bulge/bar enhancement. However, unlike Type II 
discs, inner anti-truncated discs do not appear to be structurally dis¬ 
turbed relative to Type I discs. Therefore, the additional bulge and 
outer disc light does not appear to result from restructuring of inner 
disc stars. 

The transfer of angular momentum from a bar stmcture into 
the disc would cause an increase in disc scale length outside a 
break radius (i.e. outer disc stars are redistributed to higher radii; 
|Minchev et al.|20T^ . However, this mechanism does not explain 
the intrinsic increase in disc luminosity relative to Type I discs. In¬ 
stead, it would be necessary to invoke additional star formation to 
build this additional stellar mass. If the progenitor disc was gas rich, 
outward angular momentum transfer from the bar could lead to an 
increased gas density at larger radii, and hence yield heightened 
star formation in the outer disc. Additionally, if disc gas within the 
break radius was simultaneously driven inwards by the bar, then the 
resulting central burst of star formation could explain the increased 
bulge luminosity. Central or outer starburst scenarios would be eas¬ 
ily confirmed via the optical colours of these structural components 
(i.e. systematically bluer, indicating recent star formation). How¬ 
ever, such a multi-band analysis is beyond the scope of the present 
work. 

If secular angular momentum transfer due to bar components 
is the primary mechanism of disc break formation, then the dis¬ 
tinction between Type II and III galaxies may be due to the ab¬ 
sence/presence of cold gas in the progenitor disc: a gasless (i.e. 
quenched) progenitor would result in the redistribution of disc stars, 
and hence form a Type II disc, while the bar in a gas-rich progenitor 
may interact primarily with gas, yielding a Type III disc. 

Alternatively, anti-truncated disc formation scenario via 
merger events has been proposed in |Borlaff et al.| ( |2014^ . Such 
a merger event would add mass (=luminosity) to the galaxy, and 


would grow the bulge component (oc ; Boylan-Kolchin et al. 


|2005[[van Dokkum et al.|2010[|Hilz et al.|2012^ . In this paradigm, 

the outer Type III disc corresponds to a merger remnant structure, 
while the inner disc represents the surviving progenitor disc (po¬ 
tentially stabilised by the presence of a bar). If brighter galaxies 
assembled more mass via mergers, then the observation of decreas¬ 
ing rbrk/^e,in with increasing Type III galaxy luminosity can be 
understood as a decreasing fraction of the primordial disc surviving 
increasing mass ratio mergers. 


6 SUMMARY AND CONCLUSIONS 

In this paper, we have presented detailed decomposition analy¬ 
ses (both bulge-disc and more complex, multi-component models) 
of ~630 Coma cluster galaxies (in the luminosity range —17 > 
Mg > —22) using CFHT i band imaging data. As this data is 
> 12x deeper than SDSS, fitting accuracy and reliability was sub¬ 
stantially improved relative to studies based on SDSS imaging data. 


This work focused on early-type galaxies (notably those with outer 
discs, i.e. SOs). 

The Sersic bulge -l- exponential disc decomposition analy¬ 
sis previously presented in [Head et al.| j2014| ‘Paper T), has 
been extended to a wider range of candidate models including 3- 
component and/or broken disc models. This has allowed a detailed 
re-investigation of the ~400 Coma cluster galaxies previously con¬ 
sidered to be poorly-described by an archetypal (central) bulge -l- 
(outer) disc morphology, in addition to the ~ 200 archetypal SOs in 
Paper Fs analysis sample. 

Rigorous model selection testing was implemented to ensure 
no dissonance exists between galaxy and (best-fit) model structure. 
We have investigated the structural properties beyond the simple 
bulge -I- (exponential) disk morphology, the multi-component struc¬ 
ture of classic ellipticals, and the role of galaxy bars in the evo¬ 
lution of disk-dominated galaxies. Furthermore, the properties of 
broken disk structures (Freeman Types II and III) have been con¬ 
trasted with the previously-considered (unbroken) exponential disk 
(Freeman Type I), allowing investigation of the formation mecha¬ 
nisms (and hence evolutionary history) of galaxies containing such 
structures. 

The key conclusions drawn from our analysis sample of 478 
reliably-fit Coma galaxies are as follows: 

i) 48 ± 3% of galaxies (N = 230) are well-described by a simple 
Sersic, or Sersic -l- exponential model, while 3(-l-) component mod¬ 
els are required to describe 42±3% of galaxies (N = 201). Hence, 
a wide range of complex structures are found for ETGs in Coma. 

ii) Disc breaks are detected in 26 ± 4% of archetypal (central bulge 
-I- outer disc) galaxies, with equal numbers of truncated (Freeman 
Type II; 12^.2%) and anti-truncated (Freeman Type III; 13^.2%) 
discs. This corresponds to a significantly higher truncated disc frac¬ 
tion, and lower anti-truncated disc fraction than has previously been 
detected for Virgo cluster galaxies. 

iii) Multi-component Sersic galaxies were resolved into a compact 
core (with n ~ 2), surrounded by large Gaussian-like structures. 
The total (combined) half-light radii for these multi-component 
Sersic galaxies are typically ~ 11 kpc. Thus, if these galaxies 
formed from the compact ‘red nuggets’ detected at high redshifts, 
then these objects require a factor of ~ 6x growth in size. 

iv) No significant variation in galaxy morphology or multi- 
component structure was detected with projected distance from the 
Coma cluster centre. Therefore, secular processes are responsible 
for the structural changes responsible for the formation of broken 
disc galaxies. 

v) Disc breaks are found overwhelmingly in barred galaxies (Type II: 

89lg% contain bars; Type III: contain bars), while the mi¬ 

nority of galaxies with unbroken discs also contain bars (42 ± 5%). 
In addition, broken discs (of both types) are structurally correlated 
with bar size. Galaxy bars therefore play an important role in the 
formation or stabilisation of Type II and Type III broken discs. 

vi) Type II discs may not be physically truncated. Rather, inner disc 
surface brightness may be suppressed in these structures, while the 
outer disc approximately preserves the progenitor disc properties. 
However, Type II disc luminosity trends are steeper than untrun¬ 
cated discs, suggesting luminosity truncation in fainter galaxies. 

vii) Significant growth of bulge size and luminosity implies a bulge 
enhancement origin (e.g. mergers, starbursts) for Type III galaxies, 
while the inner disc (r < rurk) remains structurally consistent with 
that of untruncated galaxies. Thus, ‘anti-truncated’ discs are likely 
to result from either radial redistribution of disc gas due to bars, or 
(disc-preserving) merger events. 
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Model selection techniques are biased by the assumption that 
the set of considered models contains the ‘true’ representation of 
the underlying data. Here, the detection of genuine broken disc 
galaxies would have been significantly distorted if only a narrow 
a range of models are considered. False positive broken disc detec¬ 
tion (i.e. the fraction of reported ‘broken disc’ galaxies revealed to 
have more complex, unbroken structures via a more detailed anal¬ 
ysis) can exceed 50% if 3-component and/or multi-Sersic models 
are not also considered. Thus, decomposition analyses require a 
sufficiently broad range of candidate models in order to ensure 
meaningful results. Accordingly, consideration of models inculd- 
ing more varied structural components (e.g. Ferrer bars, core-Sersic 
bulges; [Graham et al.|2003| l may provide additional insight into the 
galaxies analysed in the present work. However, this does not com¬ 
promise the results of this study, which has explored the diversity 
of galaxy structures via the best fits from the considered range of 
models. 
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APPENDIX A: MODEL TYPE EXAMPLES 

Figure |A1| presents an illustrative example galaxy best fit by a 
Sersic-only model (upper panels), and the corresponding (overfit) 
bulge - 1 - disc model (lower panels). Surface brightness plots (/rQ 
and model residuals (image — model in counts) are included for 
both model fits (top left), as measured from the galaxy and model 
thumbnails in wedges of elliptical annuli (angle cos“^(e^), where 
e is the eccentricity of the galaxy’s target ellipse). The i band resid¬ 
ual images (including only the central quarters) are presented in 
the bottom right corners (black border) for each model fit. In ad¬ 
dition, component residual images (i.e. the residual image after all 
model components except the target component are subtracted) are 
included along the bottom in panels bordered by their pi plot line 
colours (i.e. red and blue for Sersic and exponential components 
respectively). Here, the addition of a disc component improves the 
goodness of fit (lower but this improvement is not statisti¬ 
cally significant given the increased number of fitting parameters 
(increased BIG). 

Equivalent example plots for galaxies best fit by all other 
model types (except CD, and CDd due to small sample sizes) are 
presented in Figures pk2|A8| Each best fit model (upper panels) is 
compared to its next simplest equivalent model (in terms of num¬ 
ber of model components; lower panels). Hence, best-fit BD (Fig¬ 
ure |A2| > and BS (Figure |A3^ models are compared to (underfit) S 
models; best-fit BDd (Figure [A4l >, BDD (Figure [A5| >, BSD (Figure 
|A6[ ), and BSS (Figure |A7[ l models are compared to (underfit) BD 
models; and the best-fit BSDd model (Figure |A8^ is compared to 
an (underfit) BSD model. 
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Figure Al. An example galaxy best fit by an S model (DR8 ObjID 
1237667444048527567): Surface brightness profiles {p.i), residuals (A/i^), 
and i band thumbnails for the S model (top: Sersic = red), and the corre¬ 
sponding BD model (bottom: bulge = red, disc = blue). Small images depict 
isolated model components (border colours = fii plot), and the total resid¬ 
ual (black borders). The target ellipse is noted in black in all thumbnails, 
and ID xi (major axis) and 2D BIC values are included for both models. 


APPENDIX B: DETAILS OF FITTING 
BI ID Break Parameterisation 

A simple ID (outer) profile fitting procedure was used as a pre¬ 
liminary method of disk break detection. This was used primarily 
to produce realistic input parameter values for the 2D broken disk 
model fitting (see Section|B^, but also identifies a sample of can¬ 
didate broken disk galaxies. 

Galaxy surface brightness profiles (as measured along the ma¬ 
jor axis in 45° wedges) were fit with a simple linear or broken 
linear model (analogous to exponential or broken exponential). Fit¬ 
ting was restricted to the range 3.54” < r < rsky (where rgky 
is the radius at which the total model surface brightness is equal 
to 4.94 X the sky uncertainty, following the methodology in |Er-| 
|win et al.|2012| > to avoid contamination of the surface brightness 
profile by the bulge or low level sky background uncertainty. The 
inner limit (3.54”) comes from the radius at which the bulge con¬ 
tribution, B/T(r), of an average archetypal galaxy (as determined in 
preceding chapters) drops below 1%. The outer limit is increased 
relative to the analysis presented in previous chapters to allow the 


Figure A2. An example galaxy best fit by a BD model (DR8 ObjID 
1237667444048527520). As Figure [ATj for a BD model (top: bulge = red, 
disc = blue), and the corresponding S model (bottom: Sersic = red). 


outer regions of galaxy surface brightness profiles to be charac¬ 
terised. 

A ID BIC was used to identify cases where the additional 
degrees of freedom afforded by the profile break significantly im¬ 
proved the model goodness-of-fit. For such broken galaxies, inner 
and outer disk scale length values were calculated from the inner 
and outer slopes of the best-fit broken linear models. The break ra¬ 
dius was measured directly from the point at which the linear model 
switches from the inner to the outer slope. 

Following ID break detection, 215 galaxies (from an initial 
sample of 631 Coma cluster galaxies) were selected as candidate 
broken disks. Subsequent analysis stages also include galaxies with 
no ID-detected break, however such galaxies must use generic in¬ 
put parameter values for broken disk model fitting. 


B2 GALFIT 

B2.1 Initial Processing 

To measure the structural and photometric parameters of galaxy 
bulges and discs, galaxy decomposition has been carried out us¬ 
ing GALFIT (version 3.0.4), a 2D fitting routine ( [Peng et al.|2010[ ). 
Given a user-specified model (of arbitrary complexity), galfit 
varies parameters based on a non-linear chi-squared minimisa- 
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Figure A3. An example galaxy best fit by a BS model (SDSS DR8 ObjID 
1237667444048855252). As Figure [aT] for a BS model (top: bulge = red, 
Sersic = green), and the corresponding S model (bottom: Sersic = red). 


Figure A4. An example galaxy best fit by a BDd model (SDSS DR8 ObjID 
1237667323797504020). As Figurej^for a BDd model (top: bulge = red, 
inner/outer disc = blue/cyan), and the corresponding BD model (bottom: 
bulge = red, disc = blue). 


tion algorithm until no significant reduction in chi-squared (x,^) 
is found. The parameter values of this best-fit 2-component model 
are used to estimate the underlying structure and photometry of the 
target galaxy. 

For GALFIT’s primary data input, ~ 100" x 100" thumbnail 
images were extracted from the MegaCam image frames, centred 
on each target galaxy. Secondary data products, as derived from 
the imaging data, were used to improve fitting robustness. These 
data products are described in detail in Paper I. In brief: The lo¬ 
cal background sky and the underlying statistical noise map were 
independently-determined from each galaxy thumbnail. In addi¬ 
tion, the image point spread function (psf) was characterised from 
stars in the MegaCam fields (no further than 5' from each galaxy), 
and the zero point of the magnitude scale was calibrated using aper¬ 
ture photometry. 

Absolute rest-frame magnitudes were calculated by subtract¬ 
ing the distance modulus (m — M = 35.09), and applying galac¬ 
tic dust extinction (using [Schlafly & Finkbeineij 20njJX014 mag 
in the i band) and fc-corrections (using Chilingarian et al!]|2010[ 
IChilingarian & Zolotukhin|2012[ typically < 0.01 mag). 


B2.2 Initial Conditions 

For onr analysis the initial conditions for the multicomponent fits 
are based on the best-fit bulge -l- disc models presented in Paper I. 
The iterative build-up of model complexity from the best-fit val¬ 
ues of simpler models is the convention recommended for reliable 
results from GALFIT, and is used to provide a sensible starting 
point for the shapes (axial ratios), sizes, and intensity of additional 
model components. Hence, unlike Paper I, the model fitting pro¬ 
cedure in this work was not extended (i.e. model parameters are 
not pertnrbed and re-fit to more thoroughly investigate the parame¬ 
ter space) as such an approach becomes computationally expensive 
(and highly sensitive to parameter degeneracies) for 3+ component 
models. Thus, the results of each input model were the product of 
one GALFIT cycle and instead care was taken to generate sensible 
initial parameter values. In addition to building model complexity 
iteratively, multiple input models were generated for a single model 
type if the prior model’s components could be interpreted ambigu- 
onsly. For example, a best-fit BD model’s bnlge (or disc) structure 
can be used as the basis for the bulge, bar, or disc for an input BSD 
model. This build-up of model complexity is illustrated in Figure 

im 
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Figure A5. An example galaxy best fit by a BDD model (SDSS DR8 ObjID 
1237667443511591025). As Figurej^for a BDD model (top: bulge = red, 
disci = blue, disc2 = cyan), and the corresponding BD model (bottom: 
bulge = red, disc = blue). 


Figure A6. An example galaxy best fit by a BSD model (SDSS DR8 ObjID 
1237667323797635239). As Figurej^for a BSD model (top: bulge = red, 
bar = green, disc = blue), and the corresponding BD model (bottom: bulge 
= red, disc = blue). 


B2.3 Parameter Errors 

While GALFIT provides an estimate of the parameter errors these 
are underestimate by a large factor ( [Flaussler et al.|2007) . The for¬ 
mal calculation of the parameter errors is both complex and very 
computationally expensive, and has not been carried out in this 
study. In our analysis in Section 4 where necessary we adopt the 
approach of using the scatter about the observed trends as an upper 
estimate of the statistical uncertainties in the parameters. For ex¬ 
ample, in Figure 5a the observed scatter in the bulge Sersic index 
found at each luminosity bin is ~0.7 and hence if there is no in¬ 
trinsic scatter this is a reasonable estimate of the Sersic index error. 
In future work we will analyse mock images of galaxy with simi¬ 
lar multi-component structures found here in order to fully charac¬ 
terise the parameter uncertainties. 

B2.4 Internal Dust Attenuation 

In this paper we have not considered the possible effects of inter¬ 
nal dust attenuation on the observed photometric structures. While 
this can bias measured structural parameters, particularly at bluer 
wavebands for spiral galaxies ^Driver et al.|2007||Mdllenhoff et al.| 
|2006[ [Pastrav et al.|[20T3) , over 90% of our Coma sample are 


cluster early-type galaxies where the dust content is likely to be 
small ( [Kaviraj et al.||2012l >. The (B — R) colours of cluster red- 
sequence galaxies can be nearly fully accounted for by the observed 
spectroscopically-determined stellar population trends to within an 
rms scatter of only ~0.02 mag ( [Smith et al.|2009[ l. Such homogene¬ 
ity in colour is unlikely to occur unless the internal extinction is uni¬ 
formly small. Furthermore in our analysis highly inclined galaxies, 
and those with strong dust lanes or strong asymmetries are excluded 
(see Section [T2 ) i in order to minimise the possible effects of dust 
on our conclusions. 

B3 Surface Brightness Profile Typing 

For a multi-component system, it is convenient to describe the com¬ 
bined model in terms of its component surface brightness profiles. 
Profile types for Sersic -l- exponential models where first formalised 
in [Allen et al.j [2006^ based on which component dominates at 
r = 0, and how many times the component profiles intersect (see 
Figure [B2[ (. Type 1 profiles correspond to the archetypal central 
bulge -I- outer disc shucture of spirals and SOs. Type 2/Type 5 pro¬ 
files represent dominant discs/bulges at all radii, with sub-dominant 
bulges/discs. Conversely, the centrally-dominant Sersic component 
in Type 3 profiles re-dominates the model at large radii. These pro- 


© 2015 RAS, MNRAS 000.[T][26[ 


































































































22 J. T. C. G. Head et al. 




1237667323797504163 


1237667323797504163 


r [arcsec] 

0 5 10 15 20 25 30 35 40 45 50 


1237667323797504055 




Figure A7. An example galaxy best fit by a BSS model (SDSS DR8 ObjID 
1237667323797504163). As Figurej^for a BSS model (top: bulge = red, 
Sersicl = green, Sersic2 = blue), and the corresponding BD model (bottom: 
bulge = red, disc = blue). 


Figure A8. An example galaxy best fit by a BSDd model (SDSS DR8 ObjID 
1237667323797504055). As Figure |^for a BSDd model (top: bulge = 
red, bar = green, inner/outer disc = blue/cyan), and the corresponding BSD 
model (bottom: bulge = red, bar = green, disc = blue). 


files may be non-physical representations of more complex (3-1- 
component) systems. Profile Types 4 and 6 are equivalent to Types 
1 and 3 with the roles of the Sersic and exponential components 
swapped. As such, these inverted profiles may be symptoms of er¬ 
roneous fitting pathways, rather than true physical structures. 


B4 2D Broken Disk Model 

Fitting a broken disk structure requires a model profile with distinct 
inner and outer exponential scale radii, connected via a smooth 
transition. In galfit, this profile is implemented by linking two 
exponential disk profiles (Ein and Eout) with (hyperbolic) trunca¬ 
tion functions at some break radius. This (pixel surface brightness) 
profile can be expressed as: 


S(r) = ri(r)Ein(r) -FT 2 (r)Eout(r) (Bl) 

where Ti and T 2 are the outer and inner truncation functions avail¬ 
able for GALFIT (see |Peng et al.|20T0) . The full functional form of 
the broken disk profile is: 



BSS BSDd 

{a,b,c,d,e,f) 


a BSD d BSD 
bBSD eBSD 
cBSD fBSD 


Figure Bl. Graphical illustration of the relation between the models dur¬ 
ing multi-component decomposition. Black arrows indicate which models 
take input parameter values from the best fit of a simpler model. Blue ar¬ 
rows indicate models which also take input parameter values from external 
sources. Models with multiple input variants (differing in their interpreta¬ 
tion of progenitor model components) are noted in red. The inset illustrates 
multiple input generation for BSD models from the best fit BD components. 
ID fitting is described in Appendix|B 1 1 


© 2015 RAS, MNRAS 000,[Tp6| 















































































Beyond Sersic + exponential in Coma 23 



Figure B2. Cartoon surface brightness plots for Sersic (red) + disc (blue) 
systems of each Allen type jAllen et al.|2006t . TYP^ ^ profiles ai‘e termed 
‘archetypal’, while all other profiles ai‘e described as ‘atypical’. Profile 
Types 4 and 6 are inversions of Types 1 and 3 (respectively), and may indi¬ 
cate erroneous fitting results. 



Figure B3. Cartoon example of the broken disk profile, indicating surface 
brightnesses of the inner (red) and outer (blue) disks (and their untrun¬ 
cated forms), r is normalised to the inner disk scale length, Rs^in (black 
dotted line), and S is normalised to the inner disk surface brightness at 
Rs, in (black dashed line). The inner (r^rk ~ outer (r^rk) trun¬ 

cation radii are indicated by red and blue dotted lines. In this example, 
^s,out — 0-8i?Sjin) and = O.Sfis^in- 


E(r) ~ 2 \ ^ ~ 


{2-B) -+5 

^brk 


tanh 


(2-B) - +B 

f’brk 


“ 1 “ 1 ) SOjOutGxp 


Eo.inexp 
—r 


Rs 




s,out 


(B2) 


where Rs.in and Bs.out are the inner and outer disk scale radii, 
Eo,in and Eo,out are the (untruncated) central surface brightnesses 
of the inner and outer disks, and rbrk is the break radius. Here, 
Tbrk is defined as the radius at which the inner and outer disk 
surface brightnesses are 1% and 99% of their untruncated val¬ 
ues respectively. Dimensionless parameter B is defined as B = 
2.65 — 4.98 where Ar is the break softening radius (radial 

difference within which the truncated flux drops from 99% to 1%). 
An example of the broken disk profile is presented in Figure [B3] for 
a truncated (Type II) disk with a greatly exaggerated Ar. 

The surface brightness of this model component can be fully 


described by a single galfit input parameter: surface brightness 
at the break radius, /i(r = rbrk). The value of ^(r = rbrk) is 
constrained to be identical for the inner and outer disk structures, 
ensuring continuity of the total component profile. Additionally, 
the axis ratios and position angle parameters of both disks are cou¬ 
pled for structural consistency, and Ar is fixed at 0.1 pixel (0.02"). 
Hence, the broken disk profile includes only two more free fitting 
parameters (Bs,out, and rbrk; fc = 6) than the usual exponential 
disk model (k = 4; see Table]^. 

Fitting using a truncation function with galfit yields a com¬ 
ponent’s surface brightness at rbrk, rather than the total component 
magnitude. Integrating Equation |B2| to infinity, however, is non¬ 
trivial due to the tanh function. Instead the total broken disk profile 
luminosity can be approximated using: 

/■’’brk r°0 

Ltot = / Ein(r)27rdr-I- / Eout(r)27rdr (B3) 

•io -'’“brk 

which approximates the truncation as a step function at rbrk- The 
corresponding total profile magnitude is thus: 


mtot = rrizp - 2.51ogiQ[27rq] 

- 2.51ogio[Eo.inBj,i„7(2, 

-rts,in 

+ Eo,outB2,„„t - 7(2, (B4) 

\ Rs,OUt J 

where q is the common disk axis ratio, and 7 is the incomplete 
gamma function. 


B5 Results Filter 

Results filtering is applied to ensure that only galaxies which can 
be reliably characterised by (one of) the smooth, symmetric candi¬ 
date models are considered for analysis. Model selection (i.e. the 
identification of the most statistically meaningful candidate model) 
is also a key function of this filter. This process is based on the sam¬ 
ple filtering in Paper I, which describes a number of test parameters 
in greater detail (notably A, Ares, and /mask). 

The filtering process for multi-component fits is illustrated in 
Figure [BT] and summarised as follows: 

(i) Galaxies are excluded if contaminated by nearby sources based 
on the number of masked galaxy thumbnail pixels within the target 
ellipse (/mask(atarget, ^target) > 0.4) and within the inner quarter 
of the target ellipse (/mask(atarget/2, ^target) ^ 0.3). 

(ii) Asymmetrical galaxies are also removed (A > 0.2), but the 
threshold for removing galaxies based on 6 D-residual asymmetry 
was raised to Area > 0.5, as moderate residual asymmetry may 
simply indicate the presence of unfitted structural components. 

(iii) For models with multiple variants (e.g. S5Da-f; see FigurejBTJ, a 
single (best fit) model is selected for analysis based on a simple 
test. However, for models with broken discs (BDd, CDd, BSDd) 
model variants are excluded from consideration if 0.4rbrk > 
1.678Bs,out (i.e. the outer disc contributes less than 8% of its to¬ 
tal flux) or 0.2rbrk > 1.678Bs,in (i.e. less than 0.3% of the inner 
disc’s total flux is truncated). These cuts remove anomalous model 
structures resulting from the broken disc component fitting to un¬ 
intended structures. 

(iv) For broken disc (Dd) models, galaxies are removed if Bs.out < 
O.lrbrk as a bug in galfit’s truncation yields an additional 
(strong) central point source in this regime. 
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Figure B4. Flow chart illustrating multi-component fitting model selection for models S, BD, CD, BS, BDd, CDd, BDD, BSD, BSS, and BSDd. For profile type 
definitions, refer to Section[^and|Allen et al.|f2Q06). 
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Beyond Sersic + exponential in Coma 
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ObjID 

RA 

Dec. 

Z 


Re,tot 

Mi,i 

Re,l 

711 


PAi 

COi 

Mi.2 

Re,2 

712 

<?2 

PA2 

Mi,3 

Re,3 

713 

<?3 

PAg 

-Re, out 

rbrk 

Cl IT 

Cj/T 

CalT 

Model 

Profile 

Flag 

1237665427552927881 

194.875 

28.7 

0.024 

-18.364 

1.1 

-16.449 

0.233 

6.96 

0.8 

118.556 

0.0 

-18.16 

1.237 

1.0 

0.826 

121.312 

999.0 

999.0 

999.0 

999.0 

999.0 

999.0 

999.0 

0.171 

0.829 

0.0 

BD 

3 

0 

1237665427552927902 

195.926 

28.906 

0.022 

-19.763 

502.0 

-17.618 

0.53 

1.008 

0.955 

97.153 

0.0 

-17.089 

4.686 

0.128 

0.769 

13.057 

-19.488 

7.865 

1.0 

0.824 

79.766 

2.417 

5.914 

0.139 

0.085 

0.776 

BSDd 

312 

1 

1237665427552993436 

195.018 

28.603 

0.023 

-20.068 

3.73 

-18.354 

1.221 

1.927 

0.68 

167.424 

0.0 

-18.758 

2.917 

0.374 

0.838 

8.609 

-19.303 

6.369 

0.486 

0.706 

21.469 

999.0 

999.0 

0.206 

0.299 

0.495 

BSS 

311 

0 

1237665427552993478 

195.095 

28.574 

0.022 

-19.109 

3.3 

-18.388 

2.162 

2.124 

0.659 

29.906 

0.0 

-18.324 

4.354 

0.44 

0.858 

111.1 

999.0 

999.0 

999.0 

999.0 

999.0 

999.0 

999.0 

0.515 

0.485 

0.0 

BS 

3 

0 

1237665427553124587 

195.449 

28.66 

0.029 

-18.288 

2.143 

-18.288 

2.143 

1.649 

0.723 

137.698 

0.0 

999.0 

999.0 

999.0 

999.0 

999.0 

999.0 

999.0 

999.0 

999.0 

999.0 

999.0 

999.0 

1.0 

0.0 

0.0 

S 

5 

0 


Table Cl. The structural and photometric parameters of multi-component models fits (i band) for the entire galaxy sample. The column headings are described 
in Table|C2] This table displays the first 15 data rows only; the complete version will be made available online. 



ABlC^Ja^^ 

Figure B5. Histogram of ABICres values for well-fit, 2+ component mod¬ 
els in the filtered sample (TV = 344). This compares each best-fit model 
to the next simplest (valid) model, plotted here relative to the uncertainty 
in ABICres. oTes- A dashed green line is included to indicate the SiTres 
cutoff for model acceptance. 

(v) Additionally, broken discs with rbrk < 5” are removed, as the 
inner disc of such systems behave like point sources. 

(vi) A BIC test (see Equation[TJ is applied to select the best-fit model, 
which introduces the least extra fitting parameters. When compar¬ 
ing any two models, the least complex (lowest k) model is preferred 
unless the BIC value of the higher k model is at least ScTres lower. 
For a range of (valid) candidate models, each model is paired and 
tested (in increasing order of complexity) with all other models un¬ 
til a best fit is found. 

(vii) Models with (one or more) component-to-total ratios, C/T< 0.05 
are removed from consideration during the BIC test due to high pa¬ 
rameter uncertainty. This is similar to the B/T cut for the selection 


of Sersic-only models in Paper I, hut does not make assumptions 
regarding the preferred ‘simpler’ model. 

(viii) The Xv limit for (BIC-selected) models is lowered to Xv > 1-2, 
while galaxies are now flagged if 1.1 < < 1-2. This more criti¬ 

cal cut in model xi has been calibrated through visual examination 
of model residuals. 

(ix) Galaxies with disc/outer component axis ratios, q < 0.2 are re¬ 
moved, as multi-component decomposition cannot he meaningfully 
applied to edge-on systems. 

(x) Models with Type 4 Sersic/disc profiles (i.e. Type 4, xix, xxA) 
are removed due to swapping of the hulge/bar and disc roles of the 
structural components. 

(xi) Models with Type 4 Sersic/Sersic profiles (e.g. Type 4, 4xx) have 
their components swapped (e.g. bar and bulge swap) to maintain 
the ‘inner’ role of the bulge component (or ‘inner’/‘middle’/‘outer’ 
roles for components 1, 2, and 3 in BSS models). Galaxy models 
modified in fhis way are nof removed or flagged. 

(xii) Remaining models with 0.05 < C/T < 0.1 are flagged as unre¬ 
liable. 

B6 BIC Test Results 

The results of the BIC test used to select the most statistically mean¬ 
ingful model for a given galaxy is illustrated in Figure [B5| for all 
multi-component filtered sample galaxies (N = 344; i.e. exclud¬ 
ing asymmetric galaxies, contaminated images, and bad fits). Here, 
we plot the difference in BICres between the selected ‘best fit’ and 
the next simplest (lower k) valid model, relative to the uncertainty 
in that ABICres- A green dashed line is included to indicate the 
SCTres limit, below which a model would not be chosen over a sim¬ 
pler alternative. This plot is comparable with Figure B1 in Paper I, 
which plots ABICres for Sersic-l- disc and Sersic-only models. 

While a number of galaxy models cluster close to the selection 
limit, only ~20% of models exhibit an improvement of less than 
SCTres when compared to a less complex model. The results of the 
present work are thus insensitive to slight changes to the ABICres 
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Column name Description 


ObjID 

SDSS DR8 Object ID 

RA 

Object Right Ascension [degrees] 

Dec. 

Object Declination [degrees] 

z 

Object SDSS Redshift 


Total rest-frame magnitude 

Re,tot 

Upper limit total half-light radius [kpc] 

M,,i 

Component 1 rest-frame magnitude 

Re,l 

Component 1 half-light radius [kpc] 

ni 

Component 1 Sersic index 

91 

Component 1 axis ratio (b/a) 

PAi 

Component 1 position angle [degrees] 

COi 

Component 1 boxiness 

M,,2 

Component 2 rest-frame magnitude 

Re,2 

Component 2 half-light radius [kpc] 

n2 

Component 2 Sersic index 

92 

Component 2 axis ratio (b/a) 

PA 2 

Component 2 position angle [degrees] 

M,,3 

Component 3 rest-frame magnitude 

Re,3 

Component 3 half-light radius [kpc] 

713 

Component 3 Sersic index 

93 

Component 3 axis ratio (b/a) 

PA 3 

Component 3 position angle [degrees] 


^e,out 

fbrk 

Cl/T 

C2/T 

Ca/T 

Model 

Profile 

Flag 


Outer disk half-light radius [kpc] 
Disk break radius [kpc] 
Component 1 light fraction 
Component 2 light fraction 
Component 3 light fraction 
Best-fit model 
(B/D)[ 


) Allen et al 
Fitting flag 


2006 type 


Table |Cl[ (column descriptions in Table |C2j l. The structural param¬ 
eters of the best-fit model (indicated by ‘Model’) are presented for 
each galaxy, including values for the total luminosity and combined 
half-light radiuQ A value of 999.0 indicates a parameter is not 
present in the relevant best-fit model (e.g. disk break radius in an 
unbroken BD galaxy). Fit quality flags (‘Flag’) are explained in Ta- 
bleOl 


Table C2. This table describes the column headings for Table |Cl] present¬ 
ing multi-component fitting results in the i band. Best-fit model types are 
described in Section [^ (inner component/outer component) [Alien et al.| 
I types are described in Section]^ and fitting flags are described in 


TablelC3 


Flag code 

Description 

Condition 


0 

1 

2 

Normal fit 

Bad fit (removed) 

High chi-squared 

N/A 

See Figure B4 

IT < x2 < 


3 

Low component fraction 

Any 0.05 < C/T < 0.10 

4 

Small break radius 

J-fark < 5" 



Table C3. This table describes the multi-component fitting flag codes, as 
used in Table lCll 


selection limit. Therefore, model selection based on a BIC test is 
robust for comparing multi-component galaxy models. 

A more detailed discussion of the BIC test for model selec¬ 
tion is available in Paper I and |Head] ( |2014| l. These works provide 
further details on the formulation of Equation [T] and include com¬ 
parison of BIC-selected models with by-eye selection, and F-test 
model selection. 


APPENDIX C: FITTING RESULTS CATALOGUE 

Multi-component i band fitting results for the extended Coma clus- 7 This value is an upper bound to the true value based on the assumption 

ter sample {N = 631, including blue galaxies) are presented in that major axes of all model components are aligned on the sky. 
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